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Abstract 

We prove that the Einstein equations for a spherically symmetric spacetime 
in Standard Schwarzschild Coordinates (SSC) form a closed system of three 
ordinary differential equations for a family of self-similar expanding waves, 
and the critical {k = 0) Friedmann universe associated with the pure radiation 
phase of the Standard Model of Cosmology, is embedded as a single point in 
this family. Removing a scaling law and imposing regularity at the center, we 
prove that the family reduces to an implicitly defined one parameter family of 
distinct spacetimes determined by the value of a new acceleration parameter a, 
such that a = 1 corresponds to the Standard Model. We prove that all of the 
self-similar spacetimes in the family are distinct from the non-critical k ^ 
Friedmann spacetimes, thereby characterizing the critical A; = Friedmann 
universe as the unique spacetime lying at the intersection of these two one- 
parameter families. We then present a mathematically rigorous analysis of 
solutions near the singular point at the center, deriving the expansion of 
solutions up to fourth order in the fractional distance to the Hubble Length. 
Finally, we use these rigorous estimates to calculate the exact leading order 
quadratic and cubic corrections to the redshift vs luminosity relation for an 
observer at the center. It follows by continuity that corrections to the redshift 
vs luminosity relation observed after the radiation phase of the Big Bang can 
be accounted for, at the leading order quadratic level, by adjustment of the 
free parameter a. The third order correction is then a prediction. Since self- 
similar expanding waves represent possible time-asymptotic wave patterns for 
the conservation laws associated with the highly nonlinear radiation phase, 
we propose to further investigate the possibility that these corrections to 
the Standard Model might be the source of the anomalous acceleration of 
the galaxies, an explanation wholly within Einstein's equations with classical 
sources, and not requiring Dark Energy or the cosmological constant.^ 



*This paper fills in the proofs and extends the results quoted in the authors' PNAS article [15]. 
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1 Introduction 



The Einstein equations that describe the expansion of the Universe during the 
radiation phase of the expansion form a highly nonlinear system of coupled 
wave equations in the form of conservation laws, [8]. Such wave equations 
support the propagation of waves, and self-similar expansion waves are impor- 
tant because even when dissipative terms are neglected in conservation laws, 
the nonlinearities alone provide a mechanism whereby non-interacting self- 
similar wave patterns can emerge from general interactive solutions, via the 
process of wave interaction and shock wave dissipation, [11, 6, 7]. In this pa- 
per, which elaborates and extends the results announced in [15], we construct 
a continuous one parameter family of self-similar expanding wave solutions 
of the Einstein equations and prove that the Standard Model of Cosmology, 
during the pure radiation epoch, is embedded as a single point in the family. 
Moreover, we show that the singularity in the equations at the center of the 
spacetime can be transformed into a rest point of an autonomous system of 
ODE's, and consequently all solutions in the family come in tangent to the 
same (strongest) eigensolution. It follows that near the center, the expanding 
waves in the family look just like the critical Friedmann-Robertson- Walker 
spacetime with pure radiation sources (FRW)^, and the leading order correc- 
tions to FRW are quadratic in the fractional distance to the Hubble length. 
Our intention, then, is to explore the possibility that these corrections to 
FRW, far out from the center, could account for the anomalous acceleration 
of the galaxies without Dark Energy or the cosmological constant. We end 
the paper with a derivation of the quadratic and cubic corrections to the 
redshift vs luminosity relations that would distinguish the expanding wave 
spacetimes in the family from FRW at the end of the radiation phase of the 
Standard Model. This extends the quadratic correction recorded in [15]. We 
emphasize that these corrections to redshift vs luminosity are not due to an 
external acceleration of any kind, but rather are due soley to the displacement 
of the energy density by the wave. 

Our initial insight was the discovery of a new set of coordinates in which 
FRW, (unbarred coordinates), goes over to a standard Schwarzchild metric 

^In this paper we let FRW refer to the critical (k = 0) Friedmann-Robertson- Walker metric with equation 
of state p= |pc^, [13]. 



form, (barred coordinates), in such a way that the metric components depend 
only on the single self-similar variable f /t. From this we set out to find the 
general equations for such self-similar solutions. In this paper we prove that 
the PDE's for a spherically symmetric spacetime in Standard Schwarzchild 
coordinates (SSC) reduce, under the assumption p = |/9C^, to a system of 
three ordinary differential equations^ in the same self-similar variable f /i. 
After removing one scaling parameter and imposing regularity at the center, 
we prove that there exists implicitly within the three parameter family (of 
initial conditions), a continuous one parameter family of self-similar solutions 
of the Einstein equations that extends the FRW metric. This part, then, 
expands on and fills in the proofs of the results recorded in [15]. 

Because different solutions in the family expand at different rates, our ex- 
panding wave equations introduce an acceleration parameter a, (normalized 
so a = 1 is FRW), and suitable adjustment of this parameter will speed 
up or slow down the expansion rate. By continuity of the evolution with 
respect to parameters, it follows that suitable adjustment of the parameter 
a can account for the leading order correction associated with an arbitrary 
anomalous acceleration observed at any time after the radiation phase of 
Standard Model, c.f. [15]. The next step in our program will be to obtain 
the quadratic and cubic corrections to redshift vs luminosity induced by the 
expanding waves at present time, by evolving forward, up through the p = 
stage of the Standard Model, the corrections we derived here for the ex- 
panding wave perturbations at the end of the radiation phase. Matching the 
leading order correction to the data will fix the choice of acceleration param- 
eter, and the third order correction, at that choice of acceleration parameter, 
is then a verifiable prediction of the theory. This is a topic of the authors' 
current research. 

We first set out to look for expanding wave solutions of the Einstein equations 
assuming pure radiation sources p = \pc^^ because our starting idea was that 
decay to self-similar expansion waves would most likely have occurred back 
when the universe was filled with radiation, [17]. The idea is that the sound 
speed and modulus of genuine nonlinearity (GN) are maximal during the 



^As far as we are aware the only other nontrivial way the PDE's for metrics in Standard Schwarzschild 
Coordinates with perfect fluid sources reduce to ODE's, is the time independent case when they reduce to 
the Oppenheimer-Volkofi^ equations, [18] 



radiation phase, and by standard theory of hyperbohc conservation laws, the 
modulus of GN governs the rate of decay by shock wave dissipation, even 
when dissipative terms are neglected in the equations, (c.f. [11, 6, 7]). This 
makes the existence of a family of such self-similar solutions, given by exact 
expressions, all the more interesting. In contrast, we are not so interested 
in self-similar waves during the matter dominated phase j9 ~ 0, after the 
uncoupling of matter and radiation, because when the pressure is zero, the 
resulting equations (for dust) have a zero modulus of GN, and one should 
not expect significant decayJ 

Thus the expanding waves we found here introduce an apparent anomalous 
acceleration into the Standard Model without recourse to a cosmological con- 
stant, and we propose to further investigate the possibility that the observed 
anomalous acceleration of the galaxies might be due to the fact that we are 
looking outward into an expansion wave of some extent. This would pro- 
vide an explanation for the anomalous acceleration within classical general 
relativity without recourse to the ad hoc assumption of an unobserved Dark 
Energy, with its unphysical anti-gravitational properties. Because the ex- 
panding waves have a center of expansion when a 7^ 1, this would violate the 
so-called Copernican Principle^ at least on the scale of the expanding wave.^ 
But most importantly, we emphasize our anomalous acceleration parameter 
is not put in ad hoc, but rather is derived from first principles starting from 
a theory of non-interacting, self-similar expansion waves, waves that we have 
to believe are propagating in solutions during the radiation phase of the Big 
Bang.^ 

In this paper we give detailed proofs of the claims made in our PNAS article 
[15], and improve the redshift vs luminosity relation stated there to third or- 
der in redshift factor z. This is a significant extension, requiring in particular 



^Note that although a self-similar expanding wave created when p ~ \pc^ should evolve into a non- 
interacting expansion wave during the p « phase, there is no reason to believe that the solution would 
remain self-similar after the radiation phase. 

*The Copernican Principle, the principle that we should not lie in a special place in the universe, has 
been taken as a starting assumption in cosmology since Howard Robertson and Geoffrey Walker proved, in 
the early 1930's, that the Friedmann-Robertson- Walker spacetimes are characterized by the assumption that 
they be spatially homogeneous and isotropic about every pointy [18]. Of course, the galaxies and clusters of 
galaxies are evidence of violations of the principle on smaller scales, and so on length scales larger than the 
extent of the expanding waves, we may not have a violation of the Copernican Principle, c.f. our discussion 
in the Conclusion, and footnote [24]. 

^See [4] for a study of self-similar spacetimes in general relativity. 



the development of refined estimates near tlie center (Tlieorem 10, Section 
5), and tlie resolution of what we term the mirror problem (Section 9). In 
Section 2 we introduce the coordinate transformation from co-moving coor- 
dinates to SSC that puts FRW into self-similar form. In Section 3 we derive 
the expanding wave equations, and prove that FRW solves the equations. In 
Section 4 we derive exact canonical co-moving coordinates for the spacetimes 
in the family, and use them to prove that solutions of the expanding wave 
equations are distinct from the k ^ FRW spacetimes. In Section 5 we 
introduce a transformation of the independent variable that regularizes the 
apparent singularity at f = in the expanding wave equations, transforming 
the singularity into a rest point of an autonomous system of ODE's. Us- 
ing this we prove the existence of a one parameter family of solutions that 
come into the center along the same eigenvector as the pure FRW spacetime. 
Using this we give a mathematically rigorous analysis of the asymptotics of 
solutions in the family near the center, culminating in Theorem 10, which 
provides an expansion of the spacetime metrics up to fourth order, and the 
velocity up to fifth order, in the fractional distance to the Hubble Length. 
In Section 6 we show that (to leading orders) every spacetime in the family 
foliates into flat spacelike hypersurfaces which expand at a rate given by the 
modified scale factor R{t) = t"'. 

Interestingly, all of the spacetimes in the family have a cusp type singularity 
in the velocity at r = in SSC coordinates^^, but hke FRW, the inverse of 
the coordinate mapping that originally took FRW from co-moving to SSC 
coordinates also regularizes the cusp singularity at the center when a ^ 1. 
In these coordinates, which are only approximately co-moving when a 7^ 1, 
the spacetimes in the family can be compared with FRW, and are amenable 
to a calculation of redshift vs luminosity. In Sections 7 and 8 we use these 
coordinates together with the estimates in Section 5 to derive the correction 
to redshift vs luminosity induced by the expanding waves for an observer 
positioned at the center of the expanding wave spacetimes when a 7^ 1. Our 
starting point in Section 8 is the argument for deriving redshift vs luminos- 
ity in the case of the FRW Standard Model as outlined in [9], Section 11.8. 
Calculating the third order correction term requires solving the mirror prob- 
lem, the problem of accounting for a dimming of light from distant sources 
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Similar cusp type singularities were encountered in [2, 3]. 



due solely to the curvature of the spacetime in the expanding waves when 
a 7^ 1, an effect not present in the Standard Model, and too weak to influence 
the second order correction we quoted in [15]. The detailed analysis of the 
mirror problem is presented in Section 9. The analysis in Section 9 is based 
on valid asymptotics for the geodesic equations. This is perfectly valid, but 
is the only place in the paper where complete mathematical proofs are not 
provided. Concluding remarks are made in Section 10. 

Our final result is the following theorem, holding for the one parameter family 
of self-similar expansion waves which assume pure radiation sources, p = ^pc^. 
This extends the quadratic correction to redshift vs luminosity recorded in 
[15], to third order in redshift factor z}^ 



Theorem 1 The redshift vs luminosity relation, as measured by an observer 
positioned at the center of the expanding wave spacetimes (described by metric 
(5.165) below), is given up to third order in redshift factor z by 

^ { a^ - 1 9 (a^ - l)(3a^ + 5) o ^,^,, ^, A ,^ ^, 

di = 2ctlz^ -^^z^ + ^ -^ -z^ + 0{l)\a - 1|/ \ , (1.1) 

where di is luminosity distance, (c.f. (8.4) below), ct is invariant time since 
the Big Bang, and a is the acceleration parameter which distinguishes the 
expanding waves in the family. 



When a = 1, (1.1) reduces to the correct linear relation for the radiation 
phase of the Standard Model, [9]. The second and third terms in the bracket 
in (1.1) thus give the leading order quadratic and cubic corrections to the 
redshift vs luminosity relation when a 7^ 1, thereby improving the quadratic 
estimate (6.5) of [15]. Since the adjustable parameter (a^ — 1) appears in front 
of the leading order correction in (1.1), it follows, (by continuous dependence 
of solutions on parameters), that a quadratic correction to the redshift vs lu- 
minosity observed at a time after the radiation phase of the Standard Model, 
can be accounted for by suitable adjustment of the parameter a. The third 
order correction is then a prediction of the theory. In particular, note that 
when a > 1, the leading order corrections in (1.1) imply a blue-shifting of 



^^Note that ^/to -^ cto corrects relation (6.5) of [15] for Iq ^ 1, i.e. giving di the correct dimension of 
length. 



radiation relative to the Standard Model, as observed by astronomers in the 
supernova data, [3]. Noting the positive sign of the coefficient of the third 
order term, we observe that the third order term further increases the effect 
of the quadratic term in displacing the redshift vs luminosity relation from 
the linear relation of the Standard Model. 



2 Self-Similar Coordinates for the k = FRW Space- 
time 

We consider the Standard Model of Cosmology during the pure radiation 
phase, after inflation, modeled by a critical Friedman-Robertson- Walker met- 
ric with equation of state p = Ipc^. In this paper we refer to this metric as 
FRW. In co-moving^^ coordinates the gravitational metric tensor g takes the 
critical {k = 0) FRW form, [18], 

ds^ = -dt^ + R{tfdr^ + f^d^^ (2.1) 

where f = Rr measures arclength distance at fixed time t and R = R{t) is 
the cosmological scale factor. The Einstein equations 

G = kT, (2.2) 

for metrics of form (2.1) reduce to the system of ODE's 

H = ^p, (2.3) 

p = -3{p + p)H, (2.4) 

where G is the Einstein curvature tensor, T is the stress tensor for a perfect 
fluid, 

T = {p + p)u®u + pg, (2.5) 

H is the Hubble constant 



^^Since wc work with spherically symmetric spacetimes, we say the coordinate system is co-moving if the 
radial coordinate is constant along particle paths. 
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p is the energy density of radiation and p is the radiation pressure. During 
the pure radiation epoch the Stefan- Boltzmann radiation law imphes the 
equation of state 

P = \pc^- (2.6) 

Exact expressions for the solution are given in the following theorem, which 
is a corollary of Theorem 2, [15]: 



Theorem 2 Let (2.1) solve (2.2) with equation of state (2.6). Then, (as- 
suming an expanding universe R > 0), the solution of system (2.3), (2.4) 
satisfying R = at t = and R = 1 at t = 1 is given by. 



^P=lh (2-7) 

R{t) = v^. (2.8) 

In particular, the Hubble constant satisfies^^ 

mt) ^14. (2.9) 

Our starting point is the following theorem which gives a coordinate trans- 
formation that takes (2.1) to the SSC form, 

ds^ = -B(i, f)df + ,,- , df^ + f^dn^, (2.10) 

^ ^ A{tj) ^ ' 

such that A and B depend only on f jt. For this define the self-similarity 
variables 

f 
« = | (2.11) 

and 

C = -. (2.12) 



^■^Notc that H and r are scale independent relative to the scaling law r — > ar, R — > -R of the FRW metric 
(2.1), ci. [13]. 
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Theorem 3 Assume p = ^pc^, k = and R{t) = v^. Then the FEW metric 

ds'^ = -dt^ + R{t)^dr^ + f^dn'^, 
under the change of coordinates 

R{t)r^^ 



2t 



f = R{t)r, 

transforms to the SSC-metric 

dp df'^ 



t, 



ds' 



+ 



2jo2 



roi^-y'iO) 1-^^(0 



+ r'dn 



where 



1 u^ 



V = 



VABvP 
is the SSC velocity, which also satisfies 

c 



^0^ 



2v 



l + v 



2" 



:2.i3) 
:2.i4) 



(2.15) 



(2.16) 



:2.i7) 



;2.18) 



In particular, the Jacobian and inverse Jacobians corresponding to the map- 
ping (2.13), (2.14) <^^6 given by 



dx \ i \/t 



(2.19) 



J-'^^ = 



dx \J\ V -f A 



-'O 



;2.20) 



with 



kl = ^ov^(l-^) 



:2.2r 
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Here u = (w^, u^) denote the (f, f) components of the (unit) 4-velocity of the 
sources in SSC coordinates, and we include the constant i/jq to later account 
for the time re-scaling freedom in (2.10), c.f. (2.18), page 85 of [13]. 

Proof: Letting x = (x^, x^) = {t, r) denote FRW coordinates and x = 
{x^,x^) = (i,f) denote the transformed coordinates, use (2.13) and (2.14) to 
obtain 

(2.22) 



dt 


= ^0, 


dt 


ipor 


dt 


dr 


2 


dr 


r 


df 


= ^rt 


dt 


2^/t' 


dr 





which gives (2.19)-(2.21) upon using 

c.f. (2.14). From (2.20) we obtain the metric components Qa/s in {t,f) coor- 
dinates: 





V.o^(i-^ 



-gap==J-'gJ-'=\ ^^°Vj-^y ^ |, (2.24) 



^ 4 



which verifies (2.15) assuming (2.17). It remains then to verify (2.16) and 
(2.18). For (2.18), use the equations (2.7), (2.8) to get 



M 



and since by (2.24). 



we have 



* 'i + f 



= '00, 



V = 



y/AB 

1 u^ C 



VABu^ 2' 

which gives (2.18). To verify (2.16), note that the fluid is co- moving with 
respect to the FRW {t, r)-coordinates, which means that the 4-velocity is 
given by u = (1, 0). Thus 

12 



and thus v = -f^, as claimed in (2.16). D 



V'o2' 



We now assume p = ^pc^ and that solutions depend only on ^. In the next 
section we show how the Einstein equations for metrics taking the SSC form 
(2.10) reduce to a system of three ODE's. A subsequent lengthy calculation 
then shows that FRW is a special solution of these equations. 



3 The Expanding Wave Equations 



Putting the SSC metric ansatz (2.10) into MAPLE, (suppressing the bars), 
the Einstein equations G = kT reduce to the four partial differential equa- 
tions^^ 

= ^r^roo (3.1) 

= 'frT^' (3.2) 

= ^r^r" (3.3) 

= 2^r-^^^ (3.4) 




where 




^ _ BtAt 1 fAtY Br BA 



2A^B 2A\A J r rA 

B ( Bj. \ B Bj. Aj. 

Here we assume the stress tensor for a perfect fluid, 

T^ = {pc^ + pyu^ + pg'\ (3.5) 

where, as usual, p denotes the energy density, p the pressure, v the fluid 
velocity defined in terms of the fluid 4- velocity by (2.16), and we use the 
standard summation convention and indices are raised and lowered with the 

^^Beware that in [8], A is used for the dt^ coefficient and B for the dr^ coefficient of the metric. 
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metric, c.f. [8]. The main purpose of this section is to prove the fohowing two 
theorems. (We return to the notation of using (f, f)-coordinates to distinguish 
SSC coordinates from co-moving coordinates {t^r)): 



Theorem 4 Let ^ denote the self- similarity variable 






and let 



G = 



^ 



Vab' 



Assume thatA{^), G{^) and v{^) solve the ODE's 

4{1-A)v 






G 



1 



{3 + v^)G-4v_ 
( n-A\ 2{l + v^)G-4v 
{ \ A J {3 + v^)G-4v 



where 



{■}j, = {-2v^ ^ 2(3 - v^)vG - {3 - v'')G^} 
{■}d = {{^v^ - I) - ^vG + {3 - v^)G^} , 



and define the density by 



Kp 



{3 + v'^)G-4v f2' 



Then the metric 



ds^ = -B(ndf + -T^df^ + f^dn'^ 

MO 



(3.6) 



(3.7) 



(3.8) 
(3.9) 

(3.10) 



(3.11) 
(3.12) 



(3.13) 



(3.14) 



solves the Einstein equations (3.1)- (3. 4) with velocity v = v{^) and equation 
1 

3/ 



of state \pc^. 
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The next theorem confirms the consistency of equations (3.8)- (3. 10). 

Theorem 5 The Standard Schwarzchild coordinate form (2.15) of the Stan- 
dard Model of Cosmology during the radiation phase (FRW) is a particular 
solution of equations (3.8)-(3.13). 



Proof of Theorem 4: In [8] it was shown that on smooth solutions, 
(3.1)-(3.4) are equivalent to (3.1)-(3.3) together with DivjT^^ = 0, where 
DivjT^^ = can be written in the locally inertial form, 



M}}j+ {^^TU}_^ - -I^Ab{ITI} + a#(T»° - Tli) (3.15) 

_j_2KffrTiQ0rTin _ frTiQl\2\ _ A-rp22\ 

'^Ay-'-M-'-M K-'-M) ) ^'-'-m}-' 

(c.f. equation (4.8) of [8]). Here the Minkowski stresses T^ are defined in 
terms of the stress tensor (3.5) and the metric components A and B through 
the identities 



|tO'^{(p + ,c^)-^^,, 



ni = \I^T''' = \(p + P'^')^rz:i\ (3-16) 



ni - >" = {(p+/'c^)^+p 



and 



Tlj = T'-' = ^, (3.17) 



C.f. equation (4.6) of [8], 



The first observation to make is that when p = ap., a = const., the stresses 
T*-' as well as the Minkowski stresses T^ are all linear in p. To prove the 
theorem we must show that system (3.1)-(3.3) and (3.15) closes, and reduces 
to the system of three ODE's (3.8)-(3.10) with constraint (3.13), under the 
assumption that A, G, v and f^p are functions of the self-similarity variable 
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^. Our strategy is to show that when A{^), G{^) and v{^) are substituted into 
equations (3.1)- (3.3) and (3.15), all terms not depending on ^ can be written 
in the form f^T^, which we show are all of the form f'^p times functions 
of the velocity, thereby also becoming functions of ^ under the additional 
assumption that f^p be a function of ^. 

To carry this out, begin by substituting ^(^), G{^) and v{^) into the first 
three Einstein equations (3.1)- (3.3) and define 

S'^ = Kr^Tj^, (3.18) 



to obtain 



^A^ = il-A)-S'', (3.19) 

^A^ = -^S'\ (3.20) 

B A 



^^ = \{il-A) + S''}, (3.2i; 



where 



S'' = «^V^-^^ = «|:77^^^K3 + ^') (3.22) 



S'' = ^f'pi^^cv = ^[ ^.f^ ,, \^v (3.23) 



S'' = ^f'p-^^^^ = K [ ^f^ ^ \{l + 3v'). (3.24) 



Based on this define 



where 



S'^ = KwV'\ (3.25) 



-2 

SO that V^^ are functions of v given by 

yoo = 3 + ^2 (3.27) 

V^^ = Av (3.28) 

yii = l + 3i;2, (3.29) 

1/22 _ l-v^, (3.30) 
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c.f. (3.22)-(3.24) and (3.17). 

Solving (3.19) and (3.20) for ^^4^ and equating gives the following consistency 
condition 

G(l -A)- GS^^ = -S^\ (3.31) 

Using (3.26) and (3.7), we record the constraint (3.31) as 

(l-A)G 
Kw = .^ ^ ^.^ ^ . 3.32 

(3 + v^)G -iv ^ ^ 

Substituting (3.26) into (3.32) readily confirms that (3.32) is equivalent to 
the constraint (3.13). 

Using the constraint (3.32) we can eliminate kw from (3.25), and hence from 
equations (3. 19)- (3.21). That is, substituting (3.32) into (3.25) gives 

S'' = ,, ^^~^^^, V'', (3.33) 

so that (3.22)-(3.24) take the w and r^p independent form 



,01 _ {\-A)G 



Substituting (3.34) into equation (3.19) directly gives equation (3.8). Finally 
for equation (3.10), write 

which in hght of (3.19) and (3.21) is equivalent to 

«Q = G{l-(i^ + .^^)}. (3.37) 

Substituting S^^ = kwV^^ with the expressions for w and V^^ in (3.32) and 
(3.27)- (3.29), respectively, into (3.37), leads directly to (3.9). 
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Thus we have proven that if. A^ G and v are functions of ^, then the first 
three Einstein equations are equivalent to (3.8) and (3.9) together with the 
constraint (3.13). It remains only to prove that when A{^), G{^) and v{^) 
are substituted into equation (3.15), the relations (3.34)- (3.36) will again 
eliminate all terms not depending on ^ in such a way that the resulting 
equation reduces to (3.10). 

To start, move everything to the right hand side of (3.15) and multiply 
through by f^ and use (3.18) to obtain 

= f{S''}^^f[VABS''] -2VABS'' 

Now assume that A, G, v and f'^p (and hence S^^) are all functions of ^ = f/i, 
and then use (3.19) and (3.21) to eliminate S'^^ and S^^ in the two terms 
quadratic in S^-^ in the last line to get 



= -^^Sl^+^f.v^ABS^^\ -2^/ABS 



Expanding the derivative in the second term, collecting like terms, canceling 
the ±^^^^5^^ terms that arise, and multiplying through by ^^^ then leads 
to the expression 

= {-iSf), + {iESf},, + {^5"}^^^ (3-38) 



^2 B ] jy 



V 
) IV 

where for convenience we define 



1 Vab 

S = 5 = ^. (3-39) 



1^ 



and we include the labeled brackets for future reference. Using S'^^ = nwV^^ 
we can write the derivatives of S^ and S^ in terms of v^ and w;^, which leads 
directly to the following equation equivalent to (3.38). 

The following lemma gives — in terms of -u^, ^4, G and -u, and it remains then 
to show that substitution of this expression for ^ into (3.40) then leads to 
the equation (3.10) for v^. 



Lemma 1 It follows from (3.13) together with (3.8) that 

1 — A 






(3.42) 



Proof: Equation (3.41) follows direclty from (3.13) and (3.27)-(3.29), and 
can be written as 

(3.43) 



(3.44) 





(1 

KW = 


-A) 

D ' 




where we let 










D = V'^- 


- EV'\ 




Using this we can write 








w^ 


Dw^ 


^e 


Di 


w 


I- A 


I- A 


D 


But by (3.20) and (3.43), 










M 


4vE 





(3.45) 



I- A iD 



(3.46) 



19 



and by (3.43) and (3.39) 

Using (3.46) and (3.47) in (3.45) gives (3.42). D 

Putting (3.42) into (3.40) and replacing S'^ by their values in (3.22)- (3.24), 
we can solve the resulting equation for v^ to obtain 

{.}^8^^^ = ^^ {(^yii - yoi)Eyoi + 2Dyoi}^ 

+E^§ { (ET/ii - yoi)T/oi) + D (yoo + yii) }^ - AEDV^\ 

where 

2 {.}^ = {(^yii - y00)(2E - y) + (2 - 3Ei;)L»} . (3.48) 

Using (3.19)-(3.21) to replace ^4^ and B^ by expressions involving the un- 
knowns A^E^v we get 

2 {-Yd (o^) e^c = - (V°')' {(EV^' - VO^)E + 2D} ^ (3.49) 

+ (D + V") {(EV" - V'W') + D (y"" + V") }g - 4j^DV22. 

Now since the F^-^'s depend only on the variable f , it follows that the brackets 
{•}^, {-j^ and {-j^ in (3.49) are all quadratic polynomials in E with poly- 
nomials in V as coefficients. We now find these coefficients. For {-Yd start 
with (3.48) and obtain 

2 {.}^ = (-^T/Oi + 2y00) + (^yii _ 3i/00)£; ^ ^g^yOi _ 2yii)£;2^ 

which upon using (3.27)-(3.29) gives 

{■}l = (2(3 - v^)) - 8vE + 2(-l + 3v^)E^ 
so that 

{■}D = E'{-}n. (3.50) 

c.f. (3.12). For l-}^ use (3.44) to write 

- {y''f {■}a = - {y''Y {(EV' - V'')E + 21/00 - EV"^'}^ , 
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and using (3.27)- (3.29) gives 

- {y'')' Da = - (4^')' {2(3 + v^) - l2vE + (1 + 3v^)E^}^ . (3.51) 
Finally 

{D + yll) {-j^ ={D + yll) {{EV' - y01)yOi ^ ^(yOO ^ yll) j 

= (yoo + yii _ ^yoi) {yoo^yoo ^ yii) _ (yoi)2yooyoi^)|^ ^ 

which upon using (3.27)-(3.29) leads to 

{D + V''){-}^ = 16{l + v^){3 + v^) (3.52) 

-32z;(3 + 2v^ + v^)E + 16i;2(3 + i;^)^^ 

Adding (3.51) and (3.52), and simplifying gives 

- {V'f {■}a +{D + V") {-Ib = {-lo + {-li E + {-h E\ (3.53) 
where 

{•}o = 16(l-^2)(3-^4) (3.54) 

{•}i = -32i;(l-i;2)(3-i;2) (3.55) 

{.J2 = 321^2(1 -^;2)_ (^3_5g) 

Comparing (3.53) with (3.11) thus gives 

- {V'?{-}a + {D + V") {•}« = 16(1 - v')E'- {}„ . (3.57) 

Putting (3.57) into (3.49) and using (3.50) and (3.30) yields the following 
equation equivalent to (3.15): 

2 ({■},, J^fE^ i^'( = 4(1 - ■"')E' {•}„ - Y^D\l - v'). (3.58) 

Multiplying (3.58) through by the factor ({-j^, jy^DE) and using , 

D = E[{3 + v^)G - 4v] , 

(c.f. (3.44), (3.27) and (3.28)), now readily verifies that (3.58), and hence 
(3.15), is equivalent to (3.10), as claimed. Thus the proof of Theorem 4 is 
complete. D 

Proof of Theorem 5: We begin by recording the following lemma which 
gives the relevant quantities associated with SSC coordinate representation 
of the FRW metric in terms of the single variable v. 
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Lemina 2 Consider the FRW spacetime as represented in SSC coordinates 
in Theorem 3. Then 

A = \-v^, (3.59) 

^^i 2v 

D = l-v^, (3.61) 

Kw = ^ (3.63) 

^. - €^ (3.64) 



2(1 -i;2^ 
2v 

1 — v'' 



Kw^ = —^ :;^^e (3.65) 

(3.66) 



Proof: This is a straightforward consequence of (3.19)-(3.21) and (2.15)- 
(2.18) of Theorem 3. D 

To prove Theorem 5, we verify the Standard Model on equation (3.38), which 
we have shown is equivalent to (3.10). (We omit the straightforward proof 
that FRW satisfies (3.8) and (3.9) in SSC.) We now use (3.59)-(3.65) to 
convert each of the five terms labeled by brackets in (3.38) into expressions 
involving v and ^v^ alone. 

{•}// = ^^^^(1 + 6^' -3^')^^^ (3.68) 



'IV - ^ 1 _ ^2)3 
.2 



{■}/F = 2 ' ' v' (3.70) 



{■},. = -v{l+v') (3.7i: 
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Thus we need only verify 

{•}/ + • •• + {-}y = 0. (3.72) 

Adding (3.67)-(3.71) we obtain 

{■}, +■■■+ {-h = {-fg^ + ""'is':^:-"'' } ;^, (3.73) 

(l-i;2)3 + (l--f;2)3 






But from (3.62) and (3.64) 






y2 



and substituting this into (3.73) and multiplying through by k /^al^ gives 
{•}/ + • • • + Dv = {-4^'(3 - v^) + (1 + v^){l + 6v^ - 3v^)}^ 

Expanding {-j^ and {-j^ gives the final result 

-{.}^ = SvO-7v' + 5v' -! = {■}„ 
thereby verifying (3.72). This completes the proof of Theorem 5. D 

We conclude that the Standard Model of Cosmology during the radiation 
phase corresponds to a solution of the expanding wave equations (3.1)- (3. 3) 
and (3.13) with parameter i/jq accounting for the time-scaling freedom of the 
SSC metric (2.10). More generally, it is not difficult to see that the time- 
scaling f — > ^/'o^ preserves solutions of (3.1)-(3.3) and the constraint (3.15)). 
The next theorem states that modulo this scaling, distinct solutions of (3.1)- 
(3.3), (3.15) describe distinct spacetimes. Thus the three equations (3.1)- 
(3.3) with the one scaling law describe a two parameter family of distinct 
spacetimes, one of which is FRW. 

Theorem 6 The replacement i ^ iljoi takes A{^) , G{^) andv{^) to A{^/iI)q), 
G{^/i/jo) and v{^/i{jo), and this scaling preserves solutions of (3.1)-(3.3), 
(3.15). Moreover, this is the only scaling law in the sense that any two 
solutions of (3.1)-(3.3), (3.15) not related by the scaling i — )► i/jQt describe 
distinct spacetimes. 
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Proof: To prove the theorem it suffices to show that the only coordinate 
transformation that takes solutions of form (3.14) to solutions of form (3.14) 
are the time-scaling transformations t — > at. Now solutions of form (3.14) are 
diagonal metrics in which the radial coordinate is taken to be f determined 
by the areas of the spheres of symmetry, so the problem is to show that the 
only coordinate transformation of the form (f, f) —>■ {t, f) taking a metric 

ds^ = -B{^)dt^ + idf2 + f'^dn^, 

to a metric 

ds^ = -B{^)dP + idf2 + f'^dn^ 

are the time-scalings t = at. Now since both metrics are diagonal with no 
cross terms we must have 

f=0(f), 

for some function 0. Moreover, since both metrics use the same radial coor- 
dinate f, the transformation must meet the condition 

Differentiating this latter condition with respect to r gives 

i A'' 

Since the left hand side is independent of f and the right hand side is not, it 
must be that both sides are constant, implying that 

(f){t) = at 

for some (positive) constant, as claimed. D 

In summary, equations (3.1)-(3.3) and (3.15) admit three (initial value) pa- 
rameters and one scaling law that describe a two parameter family of distinct 
spacetimes, one of which is FRW. In Section 5 we show that by imposing 
regularity at the center there results a further reduction to a continuous one 
parameter family of expanding wave solutions, such that one value of the pa- 
rameter corresponds to FRW, the Standard Model of Cosmology with pure 
radiation sources. 
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4 Canonical Co-moving Coordinates and Comparison 
with the k ^ FRW Spacetimes 

As a consequence of Theorem 4, we know that solutions of equations (3.8)- 
(3.10) correspond to self-similar spacetimes that solve the Einstein equations 

2 

with p = y/9. The following general theorem gives a canonical form for such 
spacetime metrics in co-moving coordinates. We use this to prove that none 
of the spacetime metrics that solve (3.8)- (3. 10) agree with a /c-FRW metric 
for any k ^ 0. 

Theorem 7 Consider a general self-similar spacetime metric of the form, 

ds^ = -B(ndP + -T^df^ + f^dn'^, (4.1) 

MO 

and let v = v{^) be an arbitrary smooth velocity field, c.f. (3.14)- That is, 
assume all functions depend only on the self- similarity variable ^ = |, and 
assume 

.-1 1 

(4.2) 

where 



is a timelike vector field that has unit length as measured by (4-V- ^et 0(^) 
andi{j{^), respectively, be solutions of the ODE's 

(4.3) 



m- 


u' 1 




uWAB 


5 


u = 


^ 1 

^ at + ^ 


d 

df 



d^ ^-y/ABv 
d^j ^- v^(l - ^2; 



di ^ - y/ABv 

Then under the coordinate mapping defined by 



(4.4) 



f = mi (4.5) 

i = ijiOl (4.6) 
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the metric (4-i) transforms to the metric 

ds^ = -dt^ + RiCfdf^ + 2Q{C)didf + f^dn^, (4.7) 

where 

C = |, (4.8) 

and the velocity transforms to 

v = 0. (4.9) 

Since -0 = and R{C) and Q{C) depend only on the self-similarity variable 
C = p V'l^) c^re co-moving coordinates in which the metric defined by (4.1) 
remains self-similar. The coordinates are canonical in the sense that the 
coefficient of —dt^ is unity in (i, f) coordinates. 

Proof: Solving (4.2) for u^ gives 

u^ = {VABv)iP, (4.10) 

and using this in the condition that w is a timelike unit vector gives 



It follows that the components of u are given by 






Now since u is assumed to be a unit vector, the condition that (t, f) be 
co-moving is equivalent to 

w^ = 0, (4.12) 

and the condition that the coefficient of —dt^ is unity is equivalent to 

u^ = 1, (4.13) 
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where 



M = W — ; + W 7-, 

dt or 

is the {t,f) representation of vector u. Now by (4.5), 

| = V. | = -«'' + ^ (4.14) 

1 = ^' | = -f0' + A (4.15) 

where "prime" denotes ^. Using this, the co- moving condition (4.12) is 

Q = u^ = ^^0 + ^w^ = i-^(t)' + 0) wV #^ (4.16) 

which by (4.11) is equivalent to (4.3). Similarly, the unity condition (4.13) 
becomes 

l = u' = %u' + %' = i-^ij' + ij)u' + ct>u\ (4.17) 

ot or 

which by (4.11) is equivalent to (4.4). Since the Jacobian derivatives (4.14), 
(4.15) are all functions of C, it follows that the {t, f) coordinate representation 
of metric (4.1) must be of the form (4.7). D 

The next theorem states that the self-similar solutions defined by equations 
(3.8)-(3.10) are distinct from the k ^ Q FRW spacetimes. 



Theorem 8 Spacetime metrics defined by solutions of equations (3.8)-(3.10) 
are distinct from the k ^ Friedmann- Roberts on- Walker spacetimes. 



Proof: The A; 7^ FRW spacetimes in co-moving coordinates are given by 

ds' = -dt' + ^^^^dr' + R{tfr'dn\ (4.18) 

1 — /cr^ 

[18] . Since by Theorem 7 we know each spacetime metric defined by a solution 
of equations (3.8)-(3.10) can be mapped to a co-moving coordinate system 
(t, f) in which the metric takes the form (4.7), in order to prove the theorem 
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it suffices to prove that there is no coordinate mapping that takes (4.18) to 
a metric of form (4.7) such that it preserves the co-moving condition. So 
assume such a mapping (t, f) — > (t, r) does exist. Now co-moving means fiuid 
trajectories move along r = const, so we must have 

r = g{f), (4.19) 

for some function g{f). Since both (4.18) and (4.7) are normalized so the 
coefficient of —dt'^ is unity, it also follows that 

t = i+h{f), (4.20) 

for some function h{f). Thus, 

dr = g'{f) df , 

dt'^ = df + 2h'{f)didf + h'iffdf^, 

implying that (4.18) transforms to the metric 



ds' = 



-dp + j ^f , + (h'f \ df^ + 2h' didf + f^dQ'^, (4.2r 

[ 1 — /cr"' J 



and our assumption is that we must have 



^'^' ^{h'A=R{0^ (4.22) 



[ 1 — kr"^ 
and 

'2h'{f) = Q(C)^ (4.23) 

for some functions R{() and Q{C), functions only of ( = j. Condition (4.23) 
immediately implies 

h'{f) = a = const., 

so (4.20) imphes 

t = i+ar, (4.24) 

for some constant a. Now denote the left hand side of (4.22) by /, and set 
S{t) = R{t)'^. Then for a contradiction it suffices to show that the left hand 
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side of (4.22) cannot be a function of ( when k ^ 0. But / is a function of ( 
if and only if 

so assuming it can we get 

(1 - kg{fy)S'g' 



fi = 



(1 - kg^Y 
1 - kg{ff)S'g' 



(1 - %2)2 ' 
and we must have 

ft (1 - kg^)S'g' 



< 



U (1 - kg^) [aS'g' + Sg"] + 2kSg{g'Y ' 
which after separating leads to 

2%c/' g" f l\ S' , , 

The left hand side of (4.25) is a function of f alone, so differentiating with 
respect to i and noting |r = 2^ = 2i7, {H = H{t) the Hubble constant), we 
must have 

or upon using af + t = t^ 

implying 

H{t) = ^, (4.26) 

for some constant Cq. Now (4.18) must satisfy the FRW-equations 

H' = ^p- k, (4.27) 

and 

p = -3(p ^p)H = -4pH. (4.28) 
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Equations (4.26) and (4.27) imply 

. = ^(1 + ..). (4.29) 

and (4.26) together with (4.27) imply 

p = Cir^^\ (4.30) 

for some constants Co and Ci. But (4.29) and (4.30) are inconsistent unless 
/c = 0, (the case in which p{t) is inverse square.) It follows, then, that / on 
the left hand side of (4.22) cannot be written as a function of C when A; 7^ 0, 
and so there is no mapping of A; 7^ FRW to co-moving form (4.7), and the 
theorem is proved. D 

As a corollary of the proof we have: 

Corollary 1 In the case of the k = FRW metric (4-18), the mappings 
(4- 19), (4-20) can be taken as 

r = ln(f), 

t = i + af, (4.31) 

for any constant a, and the co-moving form (4-21) of FRW in self- similarity 
variable ( is 

ds^ = -dp + -df^ + a didf + f^dU^. (4.32) 



5 Leading Order Corrections to the Standard Model 
Induced by the Expanding Waves 

We extend and make rigorous the asymptotic estimates in [15]. To clarify 
the issues, note first that system (3.8)- (3. 10) takes the form 

/ AiO \ ( A{0 \ 

^ G{0 ] =F\ G{0 , (5.1) 
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where ^ > 0, and whose solutions we are concerned with in a neighborhood 
of 1^ = 0. Thus consider a general n x n system of ODE's of the form 

^V^ = F{V). (5.2) 

Define the change of independent variable 

^ = e^ 

and observe that (5.7) goes over to the autonomous system of ODE's 

^. = ^ = F{U), (5.3) 

where 

U{s) = y(e^), (5.4) 

and note that for Vq G R*^, hms^_oo U{s) = Vq iff hm^^o ^(0 = ^o, in which 
case Fiyo) = 0. Then we have the following elementary lemma: 

Lemina 3 Consider a general n x n system of form (5.2) in which F is a 
smooth function of V G R", and assume that V{^) is a smooth solution of 
(5.2) in ^ > satisfying V{^) — > Vq as ^ — > 0, for some Vq E R". Then 
F{Vo) = 0. 

Proof: If ^ — > 0, then s — > — oo, so if F{Vo) ^ 0, then the flow of the 
autonomous ODE (5.3) takes small neighborhoods of Vq to neighborhoods 
disjoint from Vq, contradicting "l/(^) — > Vq. □ 

Now assume V{^) — > Vq G R" as ^ — > 0, so that by the lemma ^(Vo) = 0. 
Expanding F{V) in Taylor series about Vq gives 

F{V) = dF{Vo){V - Vo) + 0{1)\V - VoW (5.5) 

so near Vq, V{^) satisfies 

^{ViO - VoY = dFiVoWiO - Vo), (5.6) 

to leading order in \V — Vq\. Setting W = W{^) = V{^) — Vq in (5.6) gives 
the leading order system 

^W^ = AW, (5.7) 
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where ^4 is a constant n x n matrix with real entries. We restrict to the case 
that A has n real and distinct eigenvalues. Then let B denote the n x n 
matrix that diagonalizes A, so that 

5-U5 = diag{Ai,...,Aj, 

where 

Ai < A2 < • • • < A„, 

are the eigenvalues of A. Setting Z = B^W G R" then reduces (5.7) to the 
diagonal system, 

^Z'iO = B~'ABW, (5.8) 

which has the eigensolutions 

UO = ^''R^. (5.9) 

where (A^, R^) are the eigenpairs of A. From this we see that the only solutions 
V{^) of (5.2) satisfying V{^) — > Vq as .^ — > correspond to the solutions Z(^) 
that lie in the span of the eigensolutions Zi{^) corresponding to Aj > 0. 

Note in particular that the solution V{^) will only have the smoothness at 
^ = allowed by the powers of the eigenvalues in (5.9). However, in the 
special case when the positive eigenvalues A^ are all positive integers, (the 
case we find below), the span of the corresponding linearized solutions (5.9) 
will all be infinitely differentiable at ^ = 0. One can show that solutions of 
the nonlinear equations can be obtained by expanding solutions in powers of 
^^% so in the case of positive integer eigenvalues below, arbitrary smoothness 
can be assumed at ^ = 0, (c.f. Lemma 12 below). 

So assume from here on that we have a solution V{^) of a nonlinear system 
(5.2) such that V{^) -^VoEWas^^O, so that ^(^0) = 0, and A = dF{Vo) 
has n real and distinct eigenvalues. From (5.3), (5.4) note that for Vq G R", 
lims^_oo U{s) = Vo iff hm^^o ^(0 = ^o, in which case ^(Vo) = 0. Then the 
solutions U{s) of (5.3) corresponding to the positive eigensolutions (5.9) of 
the linearized system (5.8) are just solutions in the unstable manifold of the 
linearization of (5.3) at Vq, namely, 

U{s) = Vo + J2 ^^e^''R^ : a^ G R I , (5.10) 
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where the sum on i is over all non-negative eigenvalues A^ of (iF(Vo). We thus 
have the following theorem: 

Theorem 9 The solutions V{^) satisfying (5.7) together with the condition 
V{^) — > Vo G R" as ^ — > are exactly the solutions U{s) = V{e^) of (5.3) in 
the unstable manifold of the hyperbolic rest point Vq, ^(Vo) = 0. 

Consider now (3.8)-(3.10), a system of form (5.7). We know by Theorem 3 
that the Standard Model satisfies (3.8)- (3. 10) with, 

where ipQ is the constant that corresponds to a change of time-scale, c.f. 
Lemma 3. Solving the last equation for -u, (assuming the case v — )• 0+ as 
^ — > 0+), and from here on letting the subscript "1" refer to the Standard 
Model, gives 

Thus (5.11) and (5.12) define the functions 74i(^), Gi(^),i;i(^) of the Standard 
Model, and so as ^ — > in the Standard Model we have 

M0\ /i-^i(0'\ /i\ 

GiiO = M H M • ^^-^3) 

MO J \ MO J V / 

We thus look for all solutions {A{^),G{0,v{^)) of (3.8)-(3.10) that tend to 
the rest point (1, 0, 0) as ^ — > 0, and these correspond to all solutions U{s) = 
V{e^) in the unstable manifold of rest point (1,0,0) of (5.3). But note that 
(3.8)-(3.10) are undetermined at {A,G,v) = (1,0,0) without knowing the 
limit of the ratio G/v in the limit ^ — > 0. 

To remedy this, define the variable 

H = G/v, (5.14) 

and consider (3.8)-(3.10) in the nonsingular variables {v,A,H), (re-ordered 
so as to diagonalize the linearized operator at Vb, c.f. (5.33) below.) Thus 
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without confusion, from here on let V refer to the 3- vector 



(5.15) 



V 

V= I A 
H 

so that we recover G from the identity 

G = Hv. 

Again let ^ = e^ and U{s) = V{e^) so that ^V^ = U{s). Then equation (3.10), 
in terms of variables V = (-u. A, i7), becomes 



v^^^v, = -[l^]\[(3+v^)H-4]v + *^''''^^'^ 



V 



2{-} 



D. 



{3 + v^)H -4 ' 



where 



{•}iv = {•}^^^ 



{■Yn = (-2^^ + 2(3 - v')H - (3 - v')H'} , 

{■}d = {(3^' - 1) - 4^^' + (3 - ^')//V} ; 
Equation (3.8) in terms of variables V = {v, A, H) goes over to 

4(1 - A) 



As = iA^ 
and equation (3.9) goes over to 



(3 + i;2)//-4 



G.s G 



^5 = -^^e = 5^5' 



V v 



where 



G.. 



V 



-H 



l-^\ 2{l + v^)H -4: 
A ) {3 + v^)H -4 



1 , 



and 



G^„,,_^fiz^^i(3 + „.)^_4 + ^('^)»" 



-u 



V" 



.2{-} 



£>. 



(3 + i;2)//-4 



(5.16) 

(5.17) 
(5.18) 
(5.19) 

(5.20) 



(5.2i: 



(5.22) 



(5.23) 
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The constraint (3.31) then becomes 



(1-A)H 

To summarize, system (5.16)-(5.23) shall be denoted 

C/, = eV^ = d A'(Q = F(U), (5.25) 

\ H'ii) ) 

where the change of variables 

^ = e^ U{s) = V{e'), (5.26) 

transforms ^4 = ^, so that in terms of U{s), (5.25) is an autonomous system 
of three ODE's equivalent to (5.16), ((5.20) and (5.21), namely 

U, = F{U). (5.27) 

We now have that as^— >Oors— >— oo, the solution 

corresponding to the Standard Model, has the regular limit, (c.f. (5.9), 
(5.13)), 

fv,{o\ ( vm \ /0\ 

Ui{s) = V,{0^ A,{i) = l-v^{0' ^ 1 UFo. (5.28) 

We now hnearize (5.27) about the nonsingular rest point Vq = (0, 1, 2). 

For this, write (5.27) as 

V \ ( f{U) \ 

A = F{U) = g{U) , (5.29) 

hJ, \h{U)J 

where f,g,h are defined by the RHS of system (5.16), (5.20) and (5.21), 
respectively. Then we obtain directly from the second equation (5.20) that 
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Neglecting terms second order in v and terms that vanish at (0, 1,2) on the 
RHS of the first equation (5.20) gives 



df{Vo 



-dlH 



1), 



(3//i;-4i;) =(1,0,0); 



(5.3i: 



and similarly from the third equation (5.21) we obtain 



dh{VQ) = -d 



= 0, 



1-A\ 2H -4: 



^ A J 3H -4 
H{-A) f2H-4\ 



1 



2(-l) 



{3H - 4) 



Vn 



A^ \3H - 4 J 

= (0,0, -(-1 + 3// -2))^^ 
= (0,0,-3). 



^{^h^>^-'^ 



Vn 



(5.32) 



Putting (5.30)- (5.32) together yields the diagonal matrix 

/ df{Vo) \ / 1 
dF{Vo) = dgiVo) = 02 
V dh{Vo) J V -3 



(5.33) 



implying that Vq is a hyperbolic rest point of system (5.27) with positive 
eigenvalues Ai = 1, A2 = 2 and negative eigenvalue A3 = —3. 

We conclude that solutions U{s) = Vq + W{s) of the linearized equations 

W, = dF{Vo) ■ W, (5.34) 

at the rest point Vq of system (5.25), such that U{s) tends to Vq as ^ — > 0, 
s — > — 00, are precisely solutions in the 2-dimensional unstable manifold Alo 
ofFo 




Mo = ( 1 I + Span 




(5.35) 



so that 



U{s] 




(5.36) 
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and 

V{0 =1 l + ^e I • (5.37) 




It thus follows from the Stable Manifold Theorem and the Hartman-Grobman 
Theorem for hyperbolic rest points, [5, 10], that solutions U{s) of the nonlin- 
ear system (5.27) that tend to the rest point Vq as s — > — oo (like the Standard 
Model), must be precisely the solutions that lie in the 2-dimensional unsta- 
ble manifold A4 of rest point Vq for the nonlinear system (5.27), where A4 
is tangent to A^o at Vq, and such that, in a neighborhood of Vq, nonlinear 
solutions in A4 are in 1 — 1 correspondence with the linear solutions (5.35) 
in A^o, [10]. In particular, it follows from (5.28) that the Standard Model 
^i(0 = Ui{s) is the particular solution in M. that corresponds to 

Thus in general we define the two constants i/jq and a^ for the linearized 
solutions as 

so that the general solution of the linearized equations (3.8)-(3.10) at rest 
point Vo is given by 



^(0 = 


2 ^ 




^(0 = 


c 
1-- 


4 


H{^) = 


2. 





(5.38) 
(5.39) 

The condition ijjq > guarantees an expanding solution [v > 0) near ^ = 0, 
and the condition and a^ > guarantees a spacetime "outside the black hole" 
in the sense that < ^4 < 1, near ^ = 0. In this paper we are interested in 
the case i/jq > 0, a^ > but when referring to a general solution in Alo we 
formally allow i/jq, a^ to be constants in R.^^ 



^^It is interesting that equations (3.8)-(3.10) admit solutions "inside the black hole" in a neighborhood of 
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We now derive refined estimates for solutions V{^) of the fully nonlinear 
system (5.25) that lie in the invariant manifold 7W. These estimates give 
the order at which a general solution V{^) in Ai diverges from the Standard 
Model Vi(^) near the rest point Vq. For this, define A^* to be the subset of 
Ai consisting of all orbits except the unique orbit manifold corresponding to 
the (weaker) eigenvalue A = 2. Thus the Standard Model hes in A^* because 
a 7^ 0, and all orbits in Al* enter the rest point tangent to the standard 
model as s — > — oo. 

The main purpose of this section is to prove the following theorem: 

Theorem 10 LetV{^) = {v{^),A(^),H{^)) be any solution of system (5.25) 
that lies in M.* , so that lim^^o^(0 = ^o = (0,1,2), and let Vi(^) = 
(fi(^), 74i(^), Gi(^)), with Hi{^) = G{^)v{^) be the solution for the Standard 
Model given in (5.13), (c.f. (5.11), (5.12)). Then there exist real constants 
i/jQ, (^ such that the following estimates hold: 

v{0 = vi{0+^^^^i^ie + O{l)\a-l\e, (5.40) 

AiO = l-a'vliO + ^-^^^\p^4e + Oil)\a-l\\e, (5.41) 
H{0 = H,{0-^-^i^ie + O{l)\a-l\e (5.42) 



: ^ + 0(l)|a-l|r^ (5.43) 

G{0 = M + ^^^^^^o'^' + 0(l)|a-l|r, (5.44) 



and 



A{0B{04 = l-^^\^^o^' + 0(l)|a-l|r, (5.45) 

where all of the right hand sides are functions of i/jq^, and for convenience 
we let 0(1)^^ incorporate the constant i/jq in the error terms that vanish on 
the Standard Model a = 1. Moreover, for each choice ofi/jQ > 0, a^ > 0, there 
exists a solution 

ViO = (V„(«),?/«(^0, A(^0) (5.46) 



^^In this paper 0(1) denotes a function bounded as ^ — > 0, and in the next section we will, without loss 
of generality, set Vo = 1- 
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in A^* such that (5.40)-(5.45) hold, and this corresponds to the two parameter 
family of exact solutions of the Einstein equations 

ds' = -ij,^Ba{M)dt' + ■ /, ^, df^ + f^dn^, (5.47) 

where 

and B{^) is defined in (5.45), so that Aa{-) and Ba{-) are independent ofi/jQ. 
As a first comment, note that (5.45) implies tliat 

VAB ^ VmW) = ^ {l - ^^^*ii'} + 0(1)1" - ll?"- (5.48) 



Now \fAB is tlie factor tliat converts invariant velocity -u, < -u < 1, over 
to coordinate velocity v} /u^ ^ where vf' and v} are the time and radial com- 
ponents of the 4- velocity of a particle moving at speed v in SSC coordinates, 

respectively. That is, 

df u^ /-—— 

dt u^ 
Thus (5.48) tells us that near the center ^ = 0, time dilation when a 7^ 1 is 
remarkably close to time dilation in the Standard Model. 

Note also that (5.40), (5.42) imply 

and 

but they are not sufficient to determine the fourth and fifth order terms 
given in (5.44) and (5.45), respectively. For example, (5.40) and (5.42) with 
G = Hv give 

x(l + ^^o^^^ + 0(l)|a-ir) 
= Gi(l + 0(l)|a-l|^4), (5.49) 

39 



where we used Hi = Hi{^) = 2 + 0(1)^2 ^^^^ ^^ = ^^(^) = !^ + 0(1)^3 to see 
the cancehng of the second order terms in the parentheses. But (5.49) does 
not give the precise form of the fifth order terms in (5.44). Note also that 
(5.45) implies the useful relation 

B = ^ {l - "^^^ii^'] + 0(1)|« - l\e. (5.50) 

(5.51) 
In a different direction, note that (5.40) together with 

t'i(?) = ^(1 + "?) = I'M + li'ie + owe (5.52) 

(c.f. (5.11)), gives the first two non-trivial terms in v{^) and A{^) as 

"(e) = ^M + ^^^'/'i!f + o(l)?^ (5.53) 

A{i) - 1 - gy + "'^^ '/"'V r + owe. (5.54) 



Also using Gi = ^/'o^ together with (5.52) gives 



(5.55) 



Fi(0 = — = 2 - -^oT + 0(l)r, (5.56) 

Vi z 



and using this together with (5.44) in (5.42) gives 

H{0 = 2 - ^^^i^ie + 0(l)^^ (5.57) 

but the 0(l)'s in (5.53)-(5.57) do not vanish when a = 1. 

Note finally that re-scaling ipQ corresponds to re-scaling time in the nonlinear 
system (3.8)-(3.10), and this corresponds to translation in the parameter s 
in (5.27), an invariance of solutions of autonomous ODE's. Changes in the 
parameter a, however, correspond to real changes in the underlying spacetime. 

To complete the picture, the following corollary of Theorem 10 gives asymp- 
totic formulas for the density, (c.f. (3.13): 
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Corollary 2 When a = 1, the density p ^ pi of the Standard Model is given 
exactly by 



r 
3 



Kpi = -^^, (5.58) 



.2\„.2 



0i(O = 2(l-^iK' (5-59) 

where vi = vi{^) is the Standard Model velocity, c.f. (5.11). When a ^ 1, 
the density satisfies the asymptotic expression 



i^p = ^, (5.60) 

m = 01 (0 - ^(1 - «')^?(0 + ^(4a^ - 1)(1 - aM{0 + 0{l)\a - l\e. 



(5.6i: 



Proof: For (5.59), note that from (3.13) 

where we use the Standard Model identity 

(1 + vl)Hi -4 = 0. 

To derive (5.61), start from (3.13) to get 

^ 3(1 +.-)(! -A) 

^^^^ 2 f 1 + (l+'^')^-4^ 



2i7 

write 



(5.63) 



and use (5.42) to estimate 



•^ + "'^^-' ^-"'.? + 0(l)|a-l|^^ (5.64) 



2H 2 
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Putting (5.63) and (5.64) into (5.62) and using (5.61) and 

{v' - vl){l - A) = 0{l)\a - l\e 
yields 

m = {m) - 1(1 + vi){A - A,)] (i + ^-^vi) + o(i)|a - i\e 

= 01 (0 - 1(1 + vl){A - A,) + MO'-^yf (5.65) 

-1(1 + vl){A - A,)^-^v! + 0{l)\a - l\e. 

where by (5.41) the three terms in (5.65) can be estimated by 

-(1 + vl){A -Ai) = A-Ai + vl{A - Ai) 

= {a^ - ly, + 3a2(l - a^)vf + {a^ - l)vt + 0{l)\a - 1\^^ 
= -(1 - a^)vf + (1 - a2)(3a2 - 1)^4 ^ q^^^^^ _ ^^6^ 

M0'-^vj = l{l-a')vt + O{l)\a-l\e, 

-1(1 + vl){A - A,)i^vl = -1(1 - a^vt + 0{l)\a - l\e. 
Substituting these into (5.65) and cohecting powers of vi yeilds (5.61). D 

Our strategy in the proof of Theorem 10 is to first prove the theorem in the 
simpler case when we assume the 0(1) |a — 1| terms in (5.40)-(5.45) are only 
0(1) as ^ — > 0. This is simpler because we do not need to assume continuity 
or smoothness of the function 0(1) as ^ — > 0. This is the setting for Lemmas 
4- 10 below. In the final lemma, Lemma 12 below, we will argue for the 
smoothness of these 0(1) terms in the hmit ^ — > 0, sufficient to bootstrap 
from 0(1) to 0(l)|a-l|. 

To start, recall that the invariant manifold Ai for the nonlinear system (5.25) 
is tangent to the invariant manifold A^o of the linearized system (5.34), at 
the rest point Vq, so we can conclude that M. is normal to the vector (0, 0, 1) 
in {v,A,H) space at Vq, c.f. (5.35). Moreover, since solutions in M. are 
associated with eigenvalues A = 1, 2, it follows that except for the unique orbit 
associated with the unstable manifold for eigenvalue A = 2, (orbit tangent to 
(0, 1,0) at Vo), all solutions in M. come into the rest point Vq, in backward 
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time, tangent to the linearized solutions of the strongest eigenvalue, (smallest 
in backward time), A = 1. Thus Al* is the set of solution trajectories of 
Us = F{U) that come into Vq tangent to (1,0,0). Our starting point in the 
proof of Theorem 10 is thus the following corollary of the Invariant Manifold 
Theorem: 



Lemma 4 Let V{^) = {v{^),A(^),H{^)) denote any solution of (5.25) in 
M.* , so that lim^^o^(0 = ^o- Then V{^) is a smooth function away from 
^ = 0, and 



H{0 = 2 + O{l)e, 
and there exists i/jq ^ such that 



v{0 



n.Mf; 



= wm, 



lim^^ow{^) = 0. 



(5.66) 



(5.67) 



where 0(1) denotes a function smooth for ^ > and bounded as ^ ^ 0, and 
w{^) satisfies 



(5.68) 



Note that the Standard Model (5.28) satisfies (5.66), (5.67) with the much 
stronger estimate w{^) = 0(1)^^. 

We now prove a number of lemmas that improve on Lemma 4. The main 
technical result we use is the following: 



Lemma 5 Let u{^) be a solution of the scalar ODE 

^u^ = /(m,<^), 
such that u{^) is smooth for ^ > 0, and bounded as ^ — )• 0. Then: 



(5.69) 



(a) If f is a smooth function such that f{u, 0) has a finite number of non- 
degenerate rest points, then u{^) must tend to a rest point of f{u, 0). That is, 
there exists a finite uq eH such that /(O, wq) = 0, and lim.^^Qu{^) = uq. 
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(b) // equation (5. 69) is of the form 

^u'iO = auiOil + giOn. (5.70) 

where n > 0, a > 0, and g is a function continuous for ^ > and bounded 
in the limit ^ — > 0, then 

u{0 = uoCil + 0{l)C) (5.71) 

as ^ — > 0, and the sign {u{^)} is constant. 

(c) // equation (5. 69) is of the form 

^u\0 = -«w(0 + /3 + giOC, (5.72) 

where n > 0, a > 0, 13 is real, and g is a function continuous for ^ > and 
bounded in the limit ^ — > 0, then 

uiO = - + 0(l)r (5.73) 

a 

05 ^ — > 0. 

Proof of (a): Note first that (a) holds with f{u,^) replaced by /(w,0) in 
(5.69), because the substitution ^ = e* would transform the problem into 

^« = /(«,0). 

a standard scalar ODE for which it is well known that solutions tend to oo 
or to rest points as s — > — oo, ^ — > 0. Since /(w, i^ continuous in ^, it 
follows that for e sufficiently small, the function /(w, e) has a finite number 
of non-degenerate rest points that are small perturbations of the rest points 
of /(w, 0). From this it follows that solutions of 

must tend to oo or a rest point in the limit s— >— oo,^— >Oas well. 

Proof of (b): Since u{^) is smooth away from ^ = 0, and bounded as ^ — > 0, 
it follows by (a) that u{^) must tend to a rest point of the right hand side of 
(5.70) as <^ ^ 0, so 

limM(n = 0. (5.74) 
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Define y{^) for ^ > by 

uiO = ?/(0r, (5.75) 

so that (5.74) becomes 

^u'io = &iOc + Moc = «?/(or(i + o(i)e), (5.76) 

which reduces to 

y- = 0{1)C~^. (5.77) 

Integratinging (5.77) from ^ to ^, < ^ < ^, gives 

Vii) = yii) e«W«"-^"), (5.78) 

from which we conclude that y{^) -^ i/q = uq is bounded and finite as ^ — > 0. 
Thus taking ^ — > in (5.78) and replacing ^ by ^ gives 

y{0 = u, e^(i)^" = ^^0 (1 + O(l)e) . (5.79) 

Putting (5.79) into (5.75) gives (5.73) as claimed. From (5.78), the sign of 
y{^) is either everywhere non-zero or identically zero, and so sign{u{^)} is 
constant. 

Proof of (c): Since u{^) is smooth away from ^ = 0, and bounded as ^ — > 0, 
it follows again by (a) that u{^) must tend to a rest point of the right hand 
side of (5.73) as ^ — > 0, so that 

limM(0 = - = uq. (5.80) 

C^o a 

Using (5.80) in (5.73) then gives 

^{u - wo)e = -(^{u - uq) + 0(1)^". (5.81) 

But the only way (5.81) can be satisfied by a function u[^) smooth for ^ > 0, 
and hounded as ^ — > 0, is if 

^(O = ^o + 0(l)r. (5.82) 

To see this, define y{^) for ^ > by 

u-u^ = yiOr'', (5.83) 
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so that (5.82) becomes 

^u - uo)^ = ^?/'(or" - «?/(or" = -«?/(or" + o(i)e (5.84) 

which reduces to 

y'iO = 0{l)C^^-\ (5.85) 

Integrating (5.85) from to ^ then gives 

y{0 = ?/(o) + o(i)r+", 

which in (5.83) gives 

u{0 = uo + y{0)r'' + 0{l)C. (5.86) 

Since u is bounded as ^ — > 0, it follows that y{0) = in (5.86), in which case 
(5.86) gives (5.73) . D 

Now define 

^„(0 = l-aV(0, (5.87) 

so that 

A,(i)^l-"-^ + 0{l)i\, (5.88) 

and let K(0 denote the approximate solution 
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VaiO = ( ^a(0 I = I 1 - ^^4^ I , (5.89) 



Ha{0 J \ l+vi{0 

SO that 

VaiO = ^i(0 (5.90) 

l-AaiO = «'{l-^i(0}, (5.91) 

HaiO = ^i(0, (5.92) 

and hence K(0 differs from the exact Standard Model solution Vi(^) only in 

the yl-component. Also note that as a function of ^, 




■}Po 

Age 

4 



+ 0(l)^^ (5.93) 



2 

but the 0(1) in (5.93) does not vanish when a = 1. 
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Lemma 6 Let V{^) = {v{^),A{^),H{^)) be any solution of (5.25) m M* . 
Then there exists a^ G R such that 

A{i) = A,{^) + 0(1)^^ (5.94) 

(Again, we allow i/jq and a^ to denote any real numbers when discussing all 
solutions in A^, but restrict to positive values when restricting to cosmological 
mODEls that perturb the Standard Model near ^ = 0.) 

Proof: To verify (5.94), consider equation (3.8): 

^A, = ^'-^> = _^(l^iL (5 95) 

^^ {3 + v^)Hv-4v {3 + v^)H-4 ^^^ 

Now since V{^) exactly solves (5.95), substituting (5.66) and (5.67) into (5.95) 
gives 

«(1 - Ak = 2(i^^^o(f)^^2) = 2(1 - '4)(1 + owe), (5.96) 

where again 0(1) is bounded as ^ — > 0. Equation (5.96) is an equation of 
form (5.71) with u = 1 — A, a = 2, n = 2,, so part (b) of Lemma 5 implies 
that 



1 - ^(0 = ^^ + 0(1)^4^ (5.97) 

where for convenience take, 

aVo 

(no sign assumed on a^). Now substituting 
we obtain 

A{o = 1 - a\{o' + o{i)e = A^io + om\ 

as claimed in (5.94). D 

The next lemma improves (5.67) for the function v{^): 
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Lemma 7 Let V{^) = {v{^),A{^),H{^)) be any solution of (5.25) m M* 
such that (5.94) of Lemma 6 holds. Then 

v{^) = v,{i) + 0{l)e. (5.98) 

Note that since vi{Cj = |^ + 0(1)^^ it follows from (5.98) that 

v{^) = '^^ + 0{l)e^ (5.99) 

but in this case 0(1) does not vanish when a = 1. 

Proof: By Lemma 6 it suffices to prove that there exists a^ G R such that 

^(0 = (^(0,^(0,^(0) with^(0 = ^a(0 + O(l)^^ such that (5.98) holds. 
For this, consider equation (5.16) in the form 

(5.100) 
and observe that (5.67) implies that 

^(0 = o{i% 

as ^ — > 0. Using this together with (5.66) and (5.94), and applying them in 
(5.17)-(5.19) gives 

H = 2 + 0(1)^2^ (5.101) 

l-A = 0{l)e. (5.102) 

{■}d = -l + 0{l)v\ (5.103) 

{■Yn = 0{l)v\ (5.104) 

Using (5.101)-(5.104) in (5.100) gives 

^v^ = -^ + 0(1)^1 (5.105) 

To estimate (5.105), define a{^) by 

v^ = am. (5.106) 
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Putting (5.106) into (5.105) gives 

a'iO = 0(1)^, (5.107) 

which integrates to 

«(0 = «(0) + O(l)e (5.108) 

Substituting (5.108) into (5.106) and setting a(0) = ^ gives 

which, since ^^ = ^i(0 + ^(1)^^ ^^ ^^11, gives (5.98) as claimed. D 
We now use Lemmas 6 and 7 to improve the estimate (5.66) for H{^) : 

Lemma 8 Let V{^) = {v{^),A{^),H{^)) be any solution of (5.25) m M* 
satisfying (5.94) of Lemma 6 and (5.98) of Lemma 7. Then 

G{0 = M + 0(l)^^ (5.109) 

H{0 = ^ + 0(l)^^ (5.110) 

VAB = ^ + 0{l)e, (5.111) 

Wo 

where the constants 0(1) vanish on the Standard Model a = 1. 

Proof: Note first that (5.110) and (5.111) foUow from (5.109) via the identi- 
ties G = Hv and y^AB = G^, and thus since (5.109) is exact, {Gi{^) = ipo^), 
on the Standard Model, it follows that all constants 0(1) in (5.109)-(5.111) 
are zero when a = 1. . 

To verify (5.109), consider equations (3.9) and (5.21) in the form 

(5.113) 
Note now that by (5.98), 

v{0' = vl{0 + O{l)e, (5.114) 
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and that on the Standard Model Vi(^) we have the identity 

2{l + vl)Hi-4 = 0, (5.115) 

with 

H, = ^, (5.116) 

c.f. (5.11), (5.28). Also, by (5.66) we have 

H = Hi + 0{l)^\ (5.117) 

and by (5.94) 

1-A = 0{e). (5.118) 

Substituting (5.114)-(5.118) into (5.112) gives 

^G^ = G{l + 0{e)}. (5.119) 

Now putting the definition G = -/= into (5.119) gives 

(VZB)^ = 0(1)^^ (5.120) 

and integrating this from to ^, (note that \^AB = -J^ — > ^ as ^ — > 0), 
gives 



which also gives 



VAB = i-{l + 0(l)a}, (5.12i: 



^i^ = 4{l + 0(l)a}, (5.122: 



and 

1 
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^0 (1 + 0(1)^4)}. (5.123) 



Now 

H-. « 



so using (5.123) in this gives 

//=M{i + o(lX4| ^5_;L24) 

Thus since 

and i; = ^^ + 0(1)^^, we conclude from (5.124) that 

G'(o = M+o(l)^^ 

as claimed in (5.109). D 

We now can give a proof of estimate (5.40), which gives the dependence of 
v{^) at the third order in ^, a refinement of (5.99): 



Lemma 9 Let V{^) = {v{^),A{^),H{^)) he any solution of (5.25) m M*. 
Then there exists a and ipo such that 

viO = ^i(0 + ^^4e + 0(1)^^ (5.125) 



Proof: Consider again equation (5.16) in the form 



e^.--f^lf[(3 + «^)^-4]. + ^(^)" 






2{.}^J\'' ' ' (3 + ^;2)//-4 



(5.126) 



Since we now have 



E = Hi + 0(1)^^ (5.127) 

V = vi + 0{l)^^ (5.128) 

v' = vl + 0{l)e, (5.129) 

(5.130) 



we can use these in (5.17)-(5.19) to get 
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= {-2 + 2(3 - v',)H, - (3 - vf)Hf} + 0(1)^ 



{■}d = 



where 



{1}^ + 0(1K', 


(5.131) 


-2+o(lx^ 


(5.132) 


{(3z;2-l)-4//^;2 + (3-z;2)//V} 




{(3^? - 1) - 4H,vl + (3 - i;?//^?} + 0(1)^4 




{-^}D+o{i)e, 


(5.133) 


-l+o(lx^ 


(5.134) 


{1}^ = {i};.(o 




{1}^ = {i}^(o 





denote {•}%, {•}jj with {vi{^), Hi{^)) substituted for {v,H), respectively. 
Since v is 0(1)^, estimates (5.127), (5.129), (5.131) and (5.133) imply that 
we can replace H by i^i, v by -ui, {-j^ by {1}^ and {-j^ by {1}^ in (5.126) 
and incur an error no greater than 0(1)^^. That is. 



^^e = 



V, 



,2{1} 



V 



D. 



+ 



[(3 + vl)H, - 4] 



2{1]d) (3 + ^9^1-4 



(5.135) 



Vi 



+ 0(1)^^ 



Note that the first two terms on the RHS of (5.135) are just f{vi,Ai,Hi) 
with V in place of vi in the last place of the first term, and A in place of Ai 
in the second term. Substituting 

^ = ^1 + (^ - ^i), 



and 



I- A I- Ai 



A 



Ai 



+ 



l-A l-Ai 



A 



Ax 



1-^1 

Ax 



+ 



'Ax - A 

AxA 
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into (5.135), gives 



\-v\ 



[(3 + vX)E^ -^\{v- vi) 

Ai-A' 



12{1}^ 



,2{l}^; (3 + i;2)//i-4, 



AiA 



(5.136) 

+o(l)^^ 



Now by (5.94), 



so 



and thus 



l-A = a\l-Ai) + 0{l)^\ 



A-A,= {l-a\l-A,)Ai)+0{l)i\ 



Ax- A 

AiA 



(l-a^)(l-^i) 

[l-a2(l-^i)]^i 
{l-a^){l-Ax 



+ o{i)e 



(l-aV)(l-^i7 



+ o{i)e 



(5.137) 



Using this in (5.136) gives 



^«j = /(t-i.^ii.ffi; 



1 -V', 



-(I -a') 



j^ [(3 + vl)H, - 4] I (v - v{) 

4{1} ' 



1 -v. 






,2{1}^; {?> + vl)H,-A 



(5.138) 



// 



+0(1)^^ 



and use the notation that 



{1}/ = {!}/(?) 
{1}„ ^ {!}„«) 

denote the two brackets {-jj and {•}// in (5.138), with the Is indicating that 
both are evaluated on the Standard Model Vi = Vi(^). Now since v — vi 
beside {l}j and v\ beside {l}jj are both 0(1)^^, it follows that the vi terms 
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in {l}j and {1}//, being 0(1)^^, can be set equal to zero incurring an error no 
greater than 0(1)^^. For the same reason, Hi = 2 + 0(1)^^ can be replaced 
by 2 in (5.138) to error 0(1)^^. The result then is that 

{1}. = j^ + owe = -I + ome 
{i}n = ^+ owe = 2+ owe, 

where we have used (5.132) and (5.134). Putting these in (5.138) gives 



^v^ = f{vi,Ai,Hi) + {v-vi)-2{l-ayf^0{l)e. (5.139) 

To solve (5.139) for v — vi, set 

viO = MO + Hm'- (5.140) 

Substituting (5.140) into (5.136), and using that vi satisfies 

^(z;i)e = /(^i,^i,//i), 
gives 

^'(0^' + Mo^' = ^m' + 2(a' - 1)^;? + o(l)^^ 

which reduces to, (using vf/^^ = ^o/8 + 0(1)^'^), 

^w\0 = -MO + ^ + 0{l)e. (5.141) 



Equation (5.141) is an equation of form (5.73) with w = w;, n = 2, a = 2, 
and 13 = ^-^, so Part (c) of Lemma 5 implies 

w{0 = ^Y^ + O{l)e. (5.142) 

Putting (5.82) together with (5.141) and (5.80) gives 



^(o = ^l(o + ^^'+o(l)^^ 
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as claimed. D 

With the third order dependence of v estabhshed in Lemma 9, we can now 
estabhsh the dependence of A and H up to sixth and fourth order in ^, as 
claimed in (5.41) and (5.42), respectively: 

Lemma 10 Let V{^) = {v{^),A{^),H{^)) be any solution of (5.25) m A^*. 
Then there exists a and i/jq such that 

Mi) = A(o + ^°'^^~°'V y+o(l);^ {5.143) 

mo = miO-^-^fV + Oil)^, (5.144) 

where 0(1) vanishes when a = 1, c.f. (5.41), (5.42) and (5.87). 

Proof: We first establish (5.144). For this, start with (5.110) and use (5.125) 
together with H = 2 + 0(1)^^ to estimate: 

H(i) = t^+o(i)e = —. — ^ ^+ome 

V, (1 + t^Vo'^a + owe) 

= ^(i-^-i^i'ie+omA+ome 

vi \ 8i;i J 

= ffi ( 1 - , /'"' . i'f,e + o{i)A + owe 
\ 8[!f+o(i)e) J 

= ffi - i^ V'oV + o(l)c^ 

as claimed. 

To establish (5.41), assume A{^) is a solution (5.20) written in the alternative 
form 

satisfying (5.94), and define the function A{^) by 

A{0 = Aa{0 + A4e. (5.146) 
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Then using (5.42) together with the v'^ = l;^ + 0(l)^^ l-A = ^^ + 0(^4), 
and Hi = 2 + 0{f), gives 

4(1 - A) 



m-^)^ = 



[(3 + vl)H, - 4] (1 - 1(1 - a^)^ie + 0{e)) 

[(3'yU] ^^^i^^--^)^o^^^o(a) 



(3 + ^;2)//i-4 



{?> + vl)Hi-^ 

(5.147) 

Putting (5.146) into (5.147) and using 
gives 

-a^e" + 4/ic'' = -2/[vor+ ^°'^^g"'''V oy+oK^),) (s.ms) 

which simphfies to 

(_A, = -2A + ^"''^^-"'K o(e). (5.149) 

Equation (5.149) is an equation of form (5.73) with w = ^4, n = 2, a = 2, 
and (3 = ° g~° , so Part (c) of Lemma 5 imphes 

l = 5%^ + 0(a. (5.150) 

Putting (5.150) together with (5.146) gives 

as claimed. D 

We now get the dependence of G, \/ AB at orders six and five in ^, ^^, as 
claimed respectively in (5.44), (5.45), these being the first orders at which G 
and AB diverge from the Standard Model. 
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Lemma 11 Let V{^) = {v{^),A{^),H{^)) be any solution of (5.25) m M* . 
Then there exists a and i/jq such that 



Gio = M + ''^\/^ 4e + oii)e , (5.151) 

^5 = ^ {l - ^^^^oY } + 0{l)e, (5.152) 



(5.153) 



where the constants 0(1) vanish on the Standard Model a = 1. 



Again note that all constants 0(1) in (5.151) and (5.111) are zero on the 
Standard Model a = 1. . 

Proof: To verify (5.151), use 

v = v,+ i^¥i„¥ + owe = .1 (i + ^-^i'le + o(i)?^) 

in 



to obtain 



ff = M f 1 _ '^^,^e- + our) , 



vi \ 4 

c.f. (5.125), (5.110). Using this together with 

v' = vi + o{i)e 

and 



V2 

in (5.112) gives, (c.f. (5.11)), 

CG, = -G { (1^) il±^lzf,,ie + 0(l)f^ + l' 
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Now using i;i = M + 0(1)^^ vj = 0(1)^2 ^nd 1 - ^ = ^ + 0{^^) in this 
we obtain 

«Q = G{ '''<^g'"'> Vfe'' + l + 0(lK'=.} (5.154) 

To integrate (5.154) make the change of variables 

G = |^, (5.155) 

so that W = \J AB. Using this in (5.154) gives 

which simplifies to 

|| = -{^^^V'o¥ + 0(l)f}. (5.156) 

Integrating (5.156) from to ^ then produces the formula 



w{i) = iy(o) {\ - ""^^^^i^ie + o{X)e^ . (5.157) 

Squaring (5.157) then gives (5.152) in light of the fact that W = v^AB, so 
VK(0) = 4-. Equation (5.151) then follows from the identity G = -rr^- □ 



Lemina 12 All of the constants 0(1) in (5.40)-(5.45) can be taken to be 
0{l)\a-l\. 

Proof: Since all of the equations (5.40)- (5.45) are exact when a = 1 with 
0(1) = 0, to verify the 0(1) is really 0(l)|a — 1|, it suffices to prove that each 
error 0(1) represents a smooth function of a and ^ all the way into ^ = 0. This 
follows so long as we can show that any solution V{^) = {v{^),A{^),H(^)) 
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of system (5.25) can be expanded uniquely in powers of ^, with coefficients 
smooth functions of a, in a neighborhood of ^ = 0. This is true essentially 
because when F is smooth, iif is a smooth function of v and ^4 in A^, and 
solutions v{^) and A{^) in M. can be expanded in powers of ^^^ and ^^^ where 
Ai = 1 and A2 = 2 are the eigenvalues of dFiVo), c.f. (5.34). Since this is 
tedious to carry out in full detail, we demonstrate with the scalar example 

in = /(^)- (5-158) 

Let A = f'{uo) where uq is a rest point, /(wq) = 0. Then write 

00 
u = Y,bnC\ (5.159) 

and assume / has the Taylor expansion 

f{u) = f{uo) + f2 ^^ (^ - ^0)" • (5.160) 

Putting (5.159) and (5.160) into (5.158) gives 

Ai:nb„r^=/M+i:^ E«"' • (5.161) 

n=l n=0 ■ \fc=0 / 

Now since the Taylor series for / converges, it follows that you can solve 
uniquely for c„ such that 

where 

SO that using /(wo) = 0, (5.161) then becomes 

n=l n=0 

Equating powers of ^^ we can solve inductively for the coefficients 5„ depend- 
ing on the initial condition 61 



A^n6,r' = E^-^' 







6. = ^-i^o.....bn-i) ^ ^5_^g2) 

nX 
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It follows by standard theorems that if A > 0, (5.162) converges for each 
60 • The point now is that the solution will be continuous at ^ = only for 
positive A, and will have all classical derivatives if A is a positive integer; 
and if / depends on a parameter a, it is clear that each 6„ will be a smooth 
function of a as well. 

From this is is not so difficult to see that since eigenvalues for solutions in M. 
are both positive integers Ai = 1 and A2 = 2, the larger eigenvalue being an 
integral multiple of the smaller one, the solution can similarly be expanded in 
powers of ^'^^ = ^ in a neighborhood of ^ = 0, with convergence by analogous 
arguments. It follows then that all coefficients in the expansion are continuous 
functions of ^ and a, up to ^ = 0, from which it follows that each 0(1) is 
really 0(1) |a — 1| given that 0(1) vanishes on the Standard Model a = 1. D 

Proof of Theorem 10: Lemmas 9 through 12 together with (5.110) estab- 
lish that every solution of (3.8)-(3.10) in A4* enters the rest point Vq accord- 
ing to the estimates (5.40)- (5.43). Conversely, the estimates (5.40)- (5.45) are 
strong enough to conclude that V{^) comes into Vq, to leading order in ^, 
like solutions of the linearized equations. Since, by the Hartman-Grobman 
Theorem, [10], solutions of the nonlinear equations in Ai* are in 1 — 1 cor- 
respondence with solutions of the linearized system (5.34) in A^o, it follows 
that for every choice of ^/'o and a^ there exists a solution of (3.8)- (3. 10) that 
comes into rest point Vq like (5.40)-(5.45) in A4*. This completes the proof 
of Theorem 10. 

Consider now the leading order corrections to the FRW metric implied by 
the two parameter family of solutions (5.46) of equations (3.1)-(3.3), (3.15). 
Expanding solutions in ^ about ^ = 0, we have shown that, modulo the 
scaling law, one eigen-solution tends to infinity as ^ — > 0, and the other two 
satisfy A{^) — > 1, B{^) — > 1, as ^ — > 0, for each value of the parameters 
i/jQ and a. Removing the singular solution, (corresponding to the eigenvalue 
A = 3 > that blows up as s — > —00, ^ — > 0), (5.41) and (5.45) imply that 
what remains is a two parameter family of SSC spacetimes satisfying 

A{0 = (l-^^^)+0{l)\a-l\e. (5.163) 
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that reduces exactly to the FRW Standard Model when a = 1. The parameter 
i/jQ corresponds to the time rescaling symmetry of the SSC equations, and the 
parameter a is a new parameter that changes the underlying spacetimes, and 
which we call the acceleration paramter. We thus have the following theorem: 



Theorem 11 The 2-parameter family of bounded solutions (5.46) of (3.1)- 
(3.3), (3.15), that extends FRW of the Standard Model in SSC coordinates to 
the spacetime metric (5.^7), is given in terms of the two parameters i/jq and 
a, up to errors of order ^^, by 

ds' = , ''*' ^ , + , "'*'', , +f-d^\ (5.165) 



where the velocity satisfies 

v{0 = <0 + O{l)\a-l\e. (5.166) 

Here i/jq is the time-scaling parameter, a = 1 corresponds to FRW, vi{^) 
denotes the SSC velocity of the Standard Model given in (2.15)-(2.18), and 
a ^ 1 introduces a new acceleration parameter which gives the leading order 
perturbation of FRW. In particular (2.16) gives 

so (5.166) implies that the velocity v is independent of the parameter a up to 
second order in ^. 

In light of (2.15), when a = 1, (5.165) reduces exactly to the FRW metric 

MO = (1-^?), (5.167) 

1 

ro (i~ vi 

Now the SSC coordinate representation of FRW depends only on H and f, 
both of which are invariant under the scaling r —>■ ar, R —>■ R/a of the FRW 
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B{0 = -^7^— 2T- (5-168) 



metric (2.1). It follows that the SSC representation of FRW is independent 
of a, and therefore independent of our choice of scale for R{t). Thus without 
loss of generality we can assume throughout that the FRW metric is scaled 
exactly so that 

R{t) = Vt, (5.169) 

c.f. (2.1) and [13]. We can also remove the time rescaling freedom by set- 
ting i/jQ = 1.^^ We conclude that to leading order, the 1-parameter family of 
expanding wave perturbations of the FRW metric is given by 

dt diT 
ds^ = — 7 ^ + -, ^ + f^d^^ 



with fourth order errors in ^, and the velocity is given to leading order by 

V = I (5.170) 

independent of a, up to third order errors in ^. 



6 A Foliation of the Expanding Wave Spacetimes into 
Flat Spacelike Hypersurfaces with Modified Scale 
Factor R{t) = T. 

To get insight into the geometry of the spacetime metric (5.165) when a 7^ 1, 
consider the extension of the FRW {t^r) coordinate transformation (2.13)- 
(2.14) to a 7^ 1 defined by 




r 



{6.1; 
ei\. {6.2; 



^^It is interesting that we needed to keep the time scale parameter ?/;o in the analysis of the equations in 
order to cast Vb as a rest point of a system of ODE's. That is, the time scale invariance of the equations 
translates into translation in s invariance in the autonomous system for U{s). Once this is done, we are free 
to set Vo = 1- 
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A straightforward caculation shows that the metric (5.165) transforms to 
(t, r)-coordinates as 

ds^ = -dt^ + tdr'^ + f'^d^^ + a(l - a)C,dUf. (6.3) 

Metric (6.3) takes the form of a A; = Friedmann-Robertson- Walker metric 
with a small correction to the scale factor, {Ra{t) = t^^^ instead of R{i) = 
t^'^\ and a corrective mixed term. In particular, the time slices t = const. 
in (6.3) are all flat space 71^, as in FRW, and the f = const slices agree with 
the FRW metric modified by scale factor Ra{t). It follows that the t = const. 
surfaces given by (6.1), (6.2), define a fohation of spacetime into flat three 
dimensional spacelike slices. Thus when a 7^ 1, (6.3) exhibits many of the 
flat space properties characteristic of the a = 1 FRW spacetime. 



7 Expanding Wave Corrections to the Standard Model 
in Approximate Comoving Coordinates 

The metric (6.3) is not co-moving with the velocity v of (5.166), even at the 
leading order, when a 7^ 1. To obtain (an approximate) co- moving frame, 
note that from (5.170), v is independent of a up to order ^^, so it follows that 
even when a 7^ 1, the inverse of the transformation (2.13), (2.14) gives, to 
leading order in ^, a co-moving coordinate system for (5.165) in which we can 
compare the Hubble constant and redshift vs luminosity relations for (5.165) 
when a 7^ 1 to the Hubble constant and redshift vs luminosity relations for 
a = 1 FRW, as measured by (2.1) and (2.9). Thus from here on, we take 
(t, r)-coordinates to be defined by the Standard Model coordinate map to 
FRW coordinates (2.13), (2.14), which corresponds to taking a = 1 in (6.1), 
(6.2). For this map we note that f = R{t)r gives f as a function of (t, r), and 
by (2.17), (2.18) it follows that 

C = M + 0(1)^3 ^s^^O. (7.1) 

Theorem 12 The inverse of the coordinate transformation (2.13), (2.14) 
maps (5.165) over to {t^r)- coordinates as 

ds^ = FaiCf {-dt'^ + tdr^} + f^dn'^, (7.2) 
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where 

F.iif = -^f^ = 1 + («' - 1)^ + 0(l)\a - 1|C^ (7.3) 

and the SSC velocity v in (5. 40) maps to the {t^r) -velocity 

w = -'^e + 0{l)\a-l\C\ (7.4) 

Note that by (5.40), (7.1) 

w = V - vi^ 0{l)\a - l\C^. (7.5) 

Proof: Using vi = C/2, the Jacobian of the transformation from SSC coor- 
dinates {i^f) to {t^r) coordinates is given in (2.19), namely. 



J = 



dx \ ^^ Vi 



Thus letting g denote the spacetime metric (5.165), in (t,r) coordinates g 
transforms to ^^ = J^gJ where 

^ V i^oV~4 Vi [ j-^, I I Vi 



1-CV4 / -1 



1 - aW4 V t J ' 
and so neglecting errors of order \a — 1\C^, the transformed metric is 



(7.6) 



This confirms (7.2) and (7.3). 

To establish (7.4), let u = {u^,u^) be the 4- velocity of the particle in SSC 
coordinates, and u = (u^, u^) the 4-velocity in (t, r)-coordinates, so that 

1 df 1 u^ 

""" Vab^^ VAB^' ^ ' ^ 
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and similarly, 

where we used (7.2). Now let v = vi + w. We show w = w + 0{l)\a — l\C^. 
For this note that the transformation law for vectors is 

(;::)- -■ ( S ) ■ 

Let 

Ro = Vtfl^-iJo^Y (7.10) 

Ri = (-^'^oV (7.11) 

denote the top and bottom rows of J^^ | J|, c.f. (2.20). Then (5.45) implies 
that in SSC coordinates, 

VZB = i-(l + 0(l)|a-ir), 
so since ^ = 0(C), (7.8) gives 



(J)^eon..(^^'»-). (r.12) 

where 

e = 0(l)|a-l|C^ (7.13) 

Using this in (7.9) we have 



Ro ( ^^'»; M Ro( *\+ ' 

rJ') 
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7.14) 



where we used that 

i?i-(^o,^i) = (7.15) 

by (2.17), true because vi is the SSC velocity of the Standard Model a = 1, 
and the fluid is co- moving with respect to the FRW metric (2.1) by Theorem 
3. Thus, putting in the errors from (5.40) of Theorem 10, we have that 

w = w + 0{l)\a - 1|C^ = v-v, + 0{l)\a - 1|C^ = "^^C- 
This completes the proof of Theorem 12. D 

Remark: Theorem 12 implies that, neglecting errors of order 0(1) |a — 1|^^ 
in (5.41) and (5.45), the SSC spacetime g corresponding to parameter values 
{i^o,a) takes the {t,r) coordinate form (7.2). It follows, then, that any cal- 
culation based on undifferentiated metric coefficients from the approximate 
metric (7.2) gives answers correct up to order 0(1) |a — 1|^^ in the original 
metric (5.47); and any calculation based first derivatives of metric coefficients 
(7.2) gives answers correct up to order 0(1) |a — 1|^^ in the original metric 
(5.47). Thus, since geodesies involve first derivatives of the metric, estimates 
based on geodesies of (7.2) give answers correct up to order 0(1) |a — 1|^^ 
in the original metric (5.47), with one important exception. Since radial 
geodesies of (7.4) can be obtained directly from (7.2) by setting dQ = and 
ds = without going to the geodesic equations, it follows that radial lightlike 
geodesies of (7.2) agree with radial hghthke geodesies of the original metric 
(5.47) up to errors 0(1) |a — 1|C^. 

The variable ( = f/t is a natural dimensionless perturbation parameter that 
has a physical interpretation in (t, r)-coordinates because, (assuming c = 1 
or t = ct), C ranges from to 1 as f ranges from zero to the horizon distance 
in FRW, (approximately the Hubble distance c/H), a measure of the furthest 
one can see from the center at time t units after the Big Bang, [18]; that is, 

Dist . _,. 

C ~ • (7.16) 

Hubble Length 

Thus expanding in ( gives an expansion in the fractional distance to the 
Hubble length, c.f. [13]. Note also that when a = 1 we obtain the FRW 
metric (2.1), where we have used R{t) = ^/i, c.f. (5.169). 
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Now for a first comparison of tlie relative expansion at a 7^ 1 to tlie expansion 
of FRW, define the Hubble constant at parameter value a, by 

where 

RaitX) = FaiOVt, 

equals the square root of the coefficient of dr^ in (7.2). Then one can easily 
show 

We conclude that the fractional change in the Hubble constant due to the 
perturbation induced by expanding waves a ^ 1 relative to the FRW of the 
Standard Model a = 1, is given by 



8 Redshift vs Luminosity Relations and the Anoma- 
lous Acceleration 



In this section we obtain the a 7^ 1 corrections to the redshift vs luminosity 
relation of FRW up to order (^, as measured by an observer positioned at 
the center C = of the expanding wave spacetimes described by the metric 
(5.165) when a ^ \}^ (Recall that C = f jt = r j^ft measures the fractional 
distance to the horizon, c.f. (7.16).) Now redshift vs luminosity depends on 
motion of the observer, but is otherwise a coordinate independent relation. 
The physically natural coordinate system in which to do the comparison 
with FRW (a = 1) would be co-moving with respect to the sources. Thus 
we restrict to the coordinates (t, r) defined by (2.13), (2.14), in which our 
one parameter family of expanding wave spacetimes are described, to leading 
order in C, by the metric (7.2). In (t, r) coordinates the spacetime is only 
co-moving up to order C^, so we will need to incorporate the errors (from 
co-moving) below to get formulas up to order C^. Note that the approximate 

^^This is of course a theoretical relation, as the pure radiation FRW spacetime is not transparent. 
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metric (7.2) as well as the exact spacetime metrics (5.47) both reduce exactly 
to the FRW metric when a = 1, c.f. (2.1), (5.169). 

For our derivation of the redshift vs luminosity relation for (7.2) we follow 
the development in Gron-Hervik [9], page 289 ff. The redshift vs luminosity 
relation calculation in [9] in the case of the Standard Model a = 1 leads to 

di = 2toz. (8.1) 

We now generalize the argument so that it applies to the spacetime metrics 
(7.2) when a 7^ 1, assuming sources moving with arbitrary velocity w. Thus 
assume radiation of frequency Ue is emitted radially by a source moving at 
velocity w relative to the co-moving observer at [te^re), and observed at a 
later time t = to at frequency uq at the center r = of the spacetime metric 
(7.2). Let Ue denote the (intermediate) frequency of the emitted radiation as 
measured by a co-moving observer fixed at position r = rg at time t = te- In 
the Standard Model a = 1 the formula (7.4) for w reduces to w; = because 
in this case the fluid is exactly co-moving in (t, r)-coordinates, implying that 
i/g = z/g when a = 1. But when a 7^ 1, (7.4) implies w ^ and we must 
account for the case Ue 7^ i^e- 

To start, let Ag denote the wavelength of the radiation emitted at (tg, f^e) and 
Ao the wavelength received at the center C = | = 0, (that is, at f = at later 
time ^o)- Define 

Energy Emitted by Source 

L = Absolute Luminosity = — (8-2) 

Time 

. . ^ Power Recieved . , 

I = Apparent Luminosity = (8-3) 

Area 



and let 



di = Luminosity Distance = I - — - 1 (8-4) 




z = Redshift Factor = 1. (8-5) 

Ae 

Using two serendipitous properties of the metric (7.2), namely, the metric 
is diagonal in co-moving coordinates, and there is no a-dependence on the 
sphere's of symmetry, it follows that the arguments in [9], Section 11.8, can 
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be modified to give the following restatement of Theorem 1, which extends 
the results of [15]. 

Theorem 13 The redshift vs luminosity relation, as measured by an observer 
positioned at the center ( = of the spacetime described by metric (5.165), 
with velocity profile (5.166), is given up to fourth order in redshift factor z 
by 

di = 2cto I z + ^^^ + ^ \ -z^ + 0{\)\a - \\z^ \ (8.6) 

where we use the fact, (c.f. (8.16) below), that z and ( are of the same order 
as C — )• 0. 

Again, note that when a = 1, (8.6) reduces to (8.1), correct for the radiation 
phase of the Standard Model, [9]. Thus (8.6) gives the leading order quadratic 
and cubic corrections to the redshift vs luminosity relation when a 7^ 1, 
thereby improving the quadratic estimate (6.5) of [15]. Since (a^ — 1) appears 
in front of the leading order correction in (8.6), it follows, (by continuous 
dependence of solutions on parameters), that the leading order part of any 
anomalous correction to the redshift vs luminosity relation of the Standard 
Model, observed at a time after the radiation phase, can be accounted for 
by suitable adjustment of the parameter a. In particular, note that when 
a > 1, the leading order corrections in (8.6) imply a blue-shifting of radiation 
relative to the Standard Model, as observed in the supernova data, [3]. 

To establish (8.6), let P denote the energy per time (power) of radiation 
received at the mirror (of a reflecting telescope) of area A, positioned at 
the coordinate center transverse to the radial direction, the radiation being 
emitted at a distant source moving at velocity w at [te^re), and received at 
t = to, r = 0. We start with the following elementary relation, (c.f. [9] page 
289): 

p^ Ajenergy) ^^ u_oAn^ (8.7) 

Here 

A{energy) 



L = 



Are 
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is the absolute luminosity, the energy per time emitted by the source, c.f. 
(8.3); the ratio of the frequencies, given by 



J^o 



^e 



= (i+^r\ 



accounts for losses of energy due to redshifting at the source, (c.f. (8.33) 
below); the ratio of proper times satisfies 

corrects proper time change at the receiver to proper time change at the 
source, (c.f. (8.31) and (8.33) below); and finally fj, is defined to be the 
fraction of the emitted radiation received at the mirror A. In the case of the 
Standard Model a = 1, equation (11.116), page 289 of [9] gives 

fA = -Tzr-:^- (^ 



47rtor 



e 



We will need the following proposition, which establishes that when a 7^ 1, 
(8.8) holds subject to the correction factor Ca given below in (8.10). 

Proposition 1 When a ^ 1, the value of fj[ for the family of spacetime 
metrics (5.47) is given by 

fA = -^^Ca, (8.9) 

where Ca is given by 

Ca = l-"^e + 0{l)\a-l\e. (8.10) 

5 



Proposition 1 solves what we call the mirror problem. That is, it gives the 
ratio Ca of an area A of light received from a distant source at a mirror 
positioned at the origin when a 7^ 1, to the corresponding area when a = 1, 
in the limit A ^ (the limit expressing that the mirror is small relative to the 
distance to the source.) The result is important for the physical interpretation 
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of the spacetimes when a ^ \ and the discussion and proof is the topic of 
Appendix 9. 

The proof of Theorem 13 rehes on the following lemma: 

Lemma 13 Assume that radiation of frequency Up is emitted by a source 
moving at velocity w at {te, re), and observed at a later time t = to at frequency 
V{) at the center r = of the spacetime metric (7.2). Then r^ is related to to 
by 

re = -^Vto, (8.11) 

1 + 1 

where (for this section) 

c = 



re 



Vt~e' 

the luminosity distance di is given by 

di = h(l + z)^_{l + ''^e) +0(l)\a-l\C 

= *o(l + ^)C (l - ^C + ^C') + 0{l)\a - 1|C^ (8.12) 

and the redshift z observed at the origin is given by 



Postponing the proof of Lemma 13, we now give the 

Proof of Theorem 13: Since redshift vs luminosity as measured at a point 
in spacetime is a coordinate independent relation, we can obtain dg as a 
function of z by substituting w as given in (7.4) into (8.13) to get C as a 
function of z, and then substituting this expression for ( into (8.12). To 
accomplish this, start from (7.3), 

^ =l-(«'-l)? + 0(l)|a-l|C^ (8.14) 



FaiC) 
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so 


(8.1c 


\) gives 
















1 + 


z= 1 + ^ 


f 
1- 


8 ^ ; 


1 


8 


iC=')+0(l)a- 


IC^ 






= 1 + 1- 


8 


^c' + (- 


2 

a — 
16 


^c^^- 


^C=')+ 0(1)0 


.-ic^ 


Thus 




















l+z=l+ 


c 

2 


2 1 
8 ^ 


— 


3, 2 

16 (" 


-l)C^ + 0{l)a 


-ic*. 



(8.15) 

In particular (8.15) implies C = 0(z), and using this we can solve (8.15) for 
C as a function of z. First 

z = ^ + 0(l)|a2-l|Cl (8.16) 

Putting this in (8.15) with C = 0(z) gives 

~' = ^-^^^ + 0{l)|a-l|C^ 
so that 

C = 2z + (a^ - 1)^2 + 0(l)|a - 1|C^ (8.17) 

Using (8.17) into (8.15) gives 



z 



^ "" ^z'-''^{a' + 2)z^ + 0{l)\a-l\z\ 



2 2 2 

so that 

C = 2z + (a^ - 1)^2 + (a^ - 1) (a^ + 2) z^ + 0{l)\a - 1|/. (8.18) 

Now by (8.12), 

d, = to(l + z)C (l - ^C + ^C) + 0{l)\a - 1|C^ (8.19) 

and using (8.18) we can express the right hand side as 

(1 + zK (l - \C + ^e'^ = (1 + ^)^ {.},{.}„ (8.20) 
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where 

{■}j = {2 + (a^ - l)z + (a^ - l)(a2 + 2)z^} , 

which upon collecting terms gives 

(1 + z> {■}! Dn = 2. {l + ^. + (f^!^if^.2| (8,21) 

+0(l)|a-l|z4 . 

Putting (8.21) into (8.19) then gives (8.6), thereby completing the proof of 
Theorem 13. D. 

It remains only to give the 

Proof of Lemina 13: We first establish equations (8.11)-(8.13). To start, 
note that the coordinate t measures geodesic time at fixed r only when a = 1, 
so define geodesic time at fixed r by 

dr = Fadt, (8.22) 

where by (7.3), Fa = Fa{() is given to leading orders by 

FaiO = ! + («'- 1)^ + \a- 1|0(0. (8.23) 

Note that the radial lightlike geodesies for metric (7.2) satisfy 

F^{-dt^ + tdr^) = 0, 
so the radial null geodesies in (t, r)-coordinates are given by 

dt ^t ^ ^ 

independent of a. Integration gives 
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re re 



and solving 

Vt~e Vto - re/2 

for Tg gives 



re = T^, (8.25) 

l + C/2 ^ ^ 



establishing (8.11' 



Consider next equation (8.12). This is a direct consequence of Proposition 1 
as follows. Equations (8.7)- (8.8) give 

^ = (l + z)Hlr%' (*^-^''' 

and using (8.26) and (8.3), the apparent luminosity £ is 

Thus the luminosity distance d£ satisfies 

where we have used the relation (8.25) between ( and rg. Substituting (8.10) 
into (8.28) and collecting terms in C then establishes (8.12). 

Finally, to estabhsh (8.13), consider a null geodesic 

dr 1 



dt v^' 

emitted from (^e, re) and received at r = a later time ^o- Then integrating 
gives 

r^ = - dr= —. (8.29) 

Now letting A^o denote the time of one period for radiation measured at a 
given frequency at t = ^o, ^ = 0, we have 

/"^° t/t _ /'*°+^^° dt 
Jte yt Jte+Me yt 
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where A^e denotes the time of one period as measured by the co-moving 
observer positioned at the emitter at time t = te- Thus 

/'^°+^^" dt _ /■*" dt _ 

Jte+Aie Vi Jte Vt 



SO 



and we therefore have 





Jte 


dt 


Ate 


Ato 





% Vto 



(8.30) 



for sufficiently small A^o- Now the frequency u associated with the period At 
measured by a co-moving observer is 

^==, (8.31) 

At 



where At is the proper time interval associated with (7.2), 



At = Fa At. 
The ratio of the emitted to received frequency at given a is then 

where we use that Fq = Fa{0) = 1, F^ = i^a(Ce) = ^aiC), and where from 
here on we use the notation 

C = Ce = ^. (8.32) 

Thus by (8.30) we conclude 



Ue Vto 



1^0 FaiOVi 
Now by (8.5), 



1 + . = ^ = ^^ = ^^^. (8.33) 
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Thus by (8.11) we have 



l + Z = 



which upon using (8.25) gives 

1 + C/2 Z/e . ^ 

We conclude that equation (8.13) follows directly from (8.34) together with 
the following lemma: 

Lemma 14 The following relation holds between the frequency v^ emitted by 
a source moving at velocity w at (tg, ^e) and the frequency V^ as measured in 
the co-moving frame at {te^re): 

^ ^ 1 + w 



Proof: Since the frequency defined in (8.31) is defined in terms of the in- 
variant time interval Ar, it is defined independent of coordinates fixed with 
the same observer. Thus to calculate v^jv^ we can assume that V^^ is the 
frequency measured in the (local) Minkowski frame fixed with an observer at 
r = ra, and u is the frequency as measured by a second observer at the same 
point but moving with velocity w with respect to the first observer. Thus the 
result derives from the change of frequency formula for Lorentz transforma- 
tions in special relativity. To derive this formula, for this argument supress 
the angular variables held constant along radial motion, and let (^, x) denote 
the (local) Minkowski frame fixed with a first observer co-moving with the 
metric at point P = [te^re), and let {t^x) be the (local) Minkowski frame 
fixed with the second observer moving with radial velocity w with respect 
to the first observer. Since Minkowski time changes agree with proper time 
changes we have 
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Now let X denote the lightlike vector displacement of one period of the radial 
lightray at (tg^^e)- Then by definition 



and 



Thus 



At = dt{X), 



At = dt{X). 



Ue dtjX) a 

V,-MiX)'a' ^^-^^^ 

where a and a give the unbarred and barred components of X, respectively, 

X = al-^—\=al- + —\ 
\dt dx ] \dt df ] 

Now the Lorentz transformation that takes barred to unbarred coordinates 
is 

0= 



Vl— w^ ^/l—w"^ 



^/l-w'^ Vl- 



'W^ 




so 



/I w \ _ 1 + w _ , 

^ = I. + /. « = /. o Q^^ (8-37) 



\l\ — w^ yl — uP- J yl 



w 



2 



and using this in (8.36) gives (8.35) as claimed, and thus the proof of (8.13), 
and hence Lemma 13, is complete. D 



9 Appendix: The Mirror Problem 



In this section we give the proof of Proposition 1, Section 8, which gives the 
ratio Ca of an area A of light received from a distance source at a mirror 
(telescope) positioned at the origin when a 7^ 1, to the corresponding area 
when a = 1, in the limit ^ — > 0, (the limit expressing that the mirror is 
small relative to the distance to the source.) To describe the mirror problem, 
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consider light emitted from a distant source located at {te, r^) and received 
at a mirror of area A positioned orthogonal to the line of sight at the center 
r = C = of our spherically symmetric expanding spacetimes a 7^ 1, at a 
later time t = to > te > 0. The problem is to determine the fraction /^ of 
the area of the 2-sphere emitting radiation at r = rg, t = ^e that reaches the 
mirror. In the case of the Standard Model a = 1, the center of the FRW 
(t, r)-coordinate system can be translated to any point. Taking the center to 
be (^e^^e), light rays leaving the source at {te^re) will follow radial geodesies 
dQ = 0. It follows that the area of the (unit) 2-sphere emitting radiation at 
r = Ve, t = te that reaches the mirror A when a = 1, is /_4 = A/4:7Ttorl, 
{r1 = t^rl when to = 1, c.f. (8.8)^^). When a 7^ 1, the 3-spaces at fixed 
time are not homogeneous and isotropic about every point like the a = 1 
FRW, and the geodesies leaving the center of a coordinate system centered 
at (^e, Te) will not follow dQ = exactly. So there is a correction factor Ca 
required in the formula (8.8) for /_4 when a ^ 1. The goal of this section is 
to prove that (8.10) gives Ca to order 0(1) |a — 1|C|, where Ce = ^e/v^-^^ 
We prove this for the approximate metric (7.2), which agrees with the exact 
spacetime metric (5.47) up to order 0(l)|a— l|Ce, c.f. Theorem 12. It suffices 
to prove the formula for Ca for the approximate metric (7.2) instead of the 
exact metric (5.47), because the argument we give below is based on the 
geodesic equations of motion which only involve first derivatives of metric 
components. That is, the correction to the redshift vs luminosity relation 
of the standard model when a 7^ 1 is given in the bracket in formula (8.6), 
and our argument will show that the first order correction to the geodesic 
equations from the standard model will yield the first order correction in the 
bracket, and the second order correction to the geodesies will yield the second 
order correction in the bracket, c.f. the discussion in the first paragraph of 

^^To sec this, consider a packet of lightlike radial geodesies covering angular area $7, emanating from an 
FRW coordinate center at r = re, and evolving up to an end at time t = to. Such curves, being radial 



lightlike geodesies (8.24), traverse the curves f = 2 (Vi — y/t^), d = Oq ^ ^, te < t < to, where f is radial 
distance measured from the new center rg , and 6 measures angles at center r^ ■ Now setting ^ = \/i — y/to, 
these curves project into the curves at time t = to given by f = 2 (^), = 9o, < ^ < re/2. Since in the 
Standard Model, t — to is flat Euclidean space, the latter curves, being at fixed time t — to, are just the 
straight lines in TZ^ emanating from center f = 0, sweeping out the angular region il at r — and the area A 
a.t f — Ve, r = 0. It thus follows that the area at the end is ^ = f^fl, where fe = R{to)re — \JtoTe is spatial 
distance at received time t = to. So the fractional area is /^ = fl/Aw = A/'^irtorl, as claimed. 

^°In this section we change notation, and set Ce equal to the fixed value C,e = fe/ \/te, the position of the 
fixed source, and we let C, ~ t j \ft denote a variable that runs from the mirror at the origin to the source, 

< C < Ce. 



Section 8). Thus in this section we neglect the C^ errors in (7.2). 

To start, note that by spherical symmetry with dQ^ = d9^ + sin^^d^^, it 
suffices to treat the single angle case of metric (7.2) with 6 = 7r/2, < < 27r, 

ds^ = FaiCf {-dt'^ + tdr'^} + tr2#^ (9.1) 

Faicr = 1 + ^c^ (9.2) 

where ( = r l^ft, and we neglect the C^ errors between (9.2) and (7.2). That 
is, let denote the angle that corresponds to the coordinate in spherical 
coordinates centered at the vertex r = r^ where the radiation is emitted, (de- 
fined precisely below). Now imagine a circular mirror of radius rg positioned 
at r = transverse to, and receiving light emitted from a star positioned at 
r = Te. Then the light that hits the circular boundary of the mirror follows 
a null geodesic that starts at r = re, ^ = i^ and ends at r = 0, t = to- 
Looking back toward the center r = from the new center at r = rg, the null 
geodesies emanating from r = r^ that hit the circular boundary of the mirror 
at r = are axially symmetric about the unique radial null geodesic (0 = 0) 
that connects t = te, r = r^ to the center of the mirror at t = to, r = 0. Thus 
by axial symmetry, the fractional change in along the specific null-geodesic 
emanating from r = Ve, t = te that hits the circular boundary of the mirror, 
specified by the condition 9 is fixed at 7r/2, will agree with the fractional 
change in angle along any of the other geodesies that leave r = r^, t = te and 
hit the circular boundary of the mirror, because they can be obtained one 
from another by rotation around the axis of symmetry joining the center of 
the mirror to the center of the star. Now by (8.9), Ca is the ratio of /_4 when 
a 7^ 1 to the ratio of fj[ when a = 1, which is just the ratio jf of the angular 
area Qq of geodesies that hit the mirror when a 7^ 1, to the angular area Qi 
of geodesies that hit the mirror when a = 1 in the limit of small angles. It 
follows from these considerations that 

Ca = ci (9.3) 

where 

Ca = lim t^, (9.4) 
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where 0o is the value of (f) at t = to, r = for a hghthke geodesic of (9.1) 
emanating from r = re, t = ^e at angle 0e, and received at the mirror at t = to, 
r = 0. (Note that a = 1 is in the numerator of (9.4) and in the denominator 
of ^ because the ratio of angular areas at r = rg of null geodesies that hit 
the same area at r = 0, is the reciprocal of the ratio of areas that get hit at 
r = by a fixed angular area at r = rg.) When a = 1, = 0e along all radial 
geodesies emanating from r = rg at angle 0e because the FRW Standard 
Model is homogeneous and isotropic about every point. So c^ does indeed 
give the ratio of the area of light received from a distant source at a mirror 
positioned at the origin when a ^ 1, to the corresponding area when a = 1, 
in the limit ^ — > 0, the definition of Ca- It remains only to estimate c^. 

Proposition 1 is a direct consequence of the following theorem, which together 
with its refinement in Theorem 15 below, is the main result of this appendix. 



Theorem 14 The constant Ca satisfies the asymptotic relation 

Ca = l-^^C! + 0{l)\a-l\Cl as fe^O, (9.5) 



where Q = r^^j y/t 



To prove Theorem 14, we must show that Ca as defined in (9.4) satisfies 
the asymptotic relation (9.5). For this, we construct the equations for the 
geodesic that starts out radially at (te,re,0e), neglect terms of order 0q, and 
expand 0o/0e in powers of Ce to obtain (9.5). Since we need to estimate 
the evolution of (\){t) along such a geodesic, it is convenient to transform 
the metric over to spherical coordinates (t, f, 0) centered at f = re — r = 0, 
assuming t = t. For this, we find, (c.f. Figure 1), 
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(r,0) 




f ^y^ 


^^^"~^-^^ r 




,/^f^ 


TT — (p 




r = a; y 


r = 


X + y = re 




Figure 1: Spherical coordinates centered at r = 


= re, f = 0. 



SO 



rsmcp = rsin0, 
rcos0 = fcos^ — Te, 



r = f — 2rer COS0 + Tg, 

f sin 
tan© = ^ . 



r cos — Tg 
Thus the coordinate transformation from (t, f, 0) to (t, r, 



r = y f^ — 2rer COS0 + Tg, 
f sine!) 



= tan 



-1 



r cos — Tf 



IS 



(9.6) 

(9.7) 



(9.8) 
(9.9) 



(9.10) 
(9.11) 

(9.12) 



Using the notation {x^,x^,x'^) = {t,r,(f)), {x^,x^,x'^) = {t,f,(j)), we compute 

dx f — Tp cos d) 

(9.13) 

(9.14) 

(9.15) 

(9.16) 

^^■" [r cos (p-Tey ^^^" ^ 

and 

dx^ -J' -^ -'^'- 
Using the elementary relation, 

sec (tan^"^ x) = — v 1 + x^, 
we get 



dx^ 


r 




dx^ 


ref sin (j) 




dx^ 


r 




dx^ 


— Te SUKJ) 


1 


dx^ 


(f COS0 — TeY 


sec^ 


dx" 


-2 _ ^^- ^Qg 


1 


dx^ 


(fcOS0 — Te)^ 


sec^ 


dx\ 


- ao t - n 1 ": 


) 



.2 

sec' 



.2 , r 



[fC0S(f) — Tg)^ 



by which we can replace (9.15), (9.16) with 

dx^ —rpsind 



dx^ ^2 

dx^ f^ — Tef cos 



dx"^ r^ 



Lemma 15 The following identities hold: 

— V + rM 

dx^ J \dx 



(9.17) 
(9.18) 



-f^V-^f^V (^-1^) 



'■■^-(£)^--^(S)^ (-) 



dx^ dx^ 2dx^ dx'^ /o or 

c^x^ c^x^ dx^ dx^ 
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Proof: Equations (9.19)-(9.21) follow directly from (9.13)-(9.18).n 

We omit the ^-component, and write the resulting 2-metric from (9.1) in 
(r, 0)-coordinates as 

g^j = diag [Flt.tr'^} , 

so that in (f, 0)-coordinates 

Lemma 16 For a, (3 = 1,2, g^p satisfies 

gn = tl {Fl - 1) I '^^^^^ I + 1 V (9.23) 



r 



922 - tUFl-lf^^^^^ + f-') (9.24) 






an = t{ (Fl - 1) (^^l£^£!^Mfi!l^ ^ . (9.25) 



Proof: By (9.22), 



2 ( /oi\2 //-a1\2 //-aO\2 



j„ = j, (§;)=* (Fl - 1) (%\ + (1^1 + r^ l'^ 



^dx^ J I '^ \dx^ J \dx^ J \dx^ J 
and using identity (9.19) gives (9.23). Similarly, 

and using identity (9.20) gives (9.24). Finally, 

^ dx^ dx^ //z?2 .dx^dx^ dx^ dx^ 2^^^^^^\ 

i/12 ga r^^1 r-1^9 iVa /""l^l n^9 O^l 0^9 O^l 0^9 I' 

ox^ ox^ I OX^ OX"^ OX^ OX"^ ox^ ox-^ J 
and using identity (9.21) gives (9.25). D 



It now follows that the metric (9.1) under transformation to rg-centered co- 
ordinates (t, f, 0) goes over to the metric g given by 



I -rl sin^ df^ + rlr^ sin^ #^ (9.26) 

+ 2(f — Ve COS 0)rer sin (f) drdcj) \ , 



where we have simplified the df"^ component using 

/- I\2 2 2 • 2 ? 

(r — Tg cos (p) = r — Tg sm 0. 

To compute Ca^, we need the equations for the null geodesic of (9.26) that 
starts out in the radial direction from the point (^e, ^e, 0e) toward the original 
origin r = where the mirror A is positioned, the point labeled (^o, ^e, 0) in 
(t, f, 0)-coordinates. To this end, using the notation [x^^x^^x^) = {t,r,(f)), 
the three geodesic equations take the form 

'^ ^ — _r'^ (Q 97A 

dX^ ~ ''dXdX ^ ^ 

d^f ^1 dx^ dx^ 

dX^ ^ dX dX ^ ^ 

d^d) ^9 dx^ dx^ 



where^^ 



dX^ '^ dX dX ' 



^'J = \9'''{-9^3,<r + 9.^,+h^,^]^ (9-30) 



are the Christoffel symbols. We restrict to null geodesies that satisfy the 
null condition ds^ = 0, which holds identically on solutions of (9.28)- (9.29) 
so long as it holds initially, [18] . We call the baseline null geodesic the radial 
null geodesic starting at t = ^e, ^ = 0, on which = 0, = 0. The condition 
d^s = when = imposes a condition independent of a, 

= -dt^ + tdr\ (9.31) 



^^Wc use the Einstein summation convention where by repeated up down indices are assumed summed, 
and indices are raised and lowered with the metric, [18]. 
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which when integrated, assuming initial conditions t = tg when f = 0, gives 
the a-independent relation 

t = U{r)= Q^ + v^] , (9.32) 

which holds on the baseline null geodesic. Since t = to when f = re, we also 
have 

^0 = Qre + VTe) , (9.33) 

for every a. 

Now although t as a function of f along the baseline null geodesies is indepen- 
dent of a, complications arise because the relation between f and the natural 
geodesic parameter A depends on a. To be precise, define the baseline null 
geodesic to be the solution 

{t,f,^) = {ta{X),ra{X),0), (9.34) 

of (9.27)-(9.29), 

taiX) = UifiX)), (9.35) 

satisfying the initial conditions 

tM = te, t(0) = SSLo = v^, (9.36) 

f„(0) = 0, ^(0) = 1, (9.37) 

c^(0) = 0, i(0) = 0. (9.38) 

Lemma 17 The baseline null geodesic fa satisfies the following identity used 
in the analysis of the (f)-geodesic below: 

f dfa\ d dfa 1 f. {o? 



1 \ (a^-l)(re-f) | 



dX J df dX y^F^ 

^ {l + 0(l)|a-l|Ce}. (9.39) 



More precisely, in the case a = 1, the baseline null geodesic is given by 



3x2/3 . -^ 



ti(A) = t J 1 + ^A 1 = ^e { 1 + ;^A + 0{l)Ct] , (9.40) 
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1/3 



1 



fi(A) = 2v^ M 1 + ^AJ - 1 ^ = A ^1 - ^-^A + 0{l)C ] , (9.41; 
and satisfies 



When a ^ 1 we have 

^(,).|l + 0(l)|«-l|Cf, (9.43) 

ra(A) = fi(A) + 0(l)|a-l|Ce'. (9.44) 

Proof: The baseline null geodesic ^^(A) satisfies the geodesic equation 

where for this argument we let Hi- denote the Christoffel symbols for the 
2-metric h taken to be 



ds'^ = -Fldt^ + Fltdr'^ = h,jdx'dx-\ 



where again 



F! = l + '\ '^'' /^' = 1 + 0(l)|a - lie, 
2 ^ {re - rf r^ = ^ 

^ t ' ^' te' 

Thus 
^1 _ 1.11 r .+.+.1-11 a(F„^)_ (a^-l)(r.-f) 

H' - ^h'U h ^h ^h 1-1 1 ^1^- 1 

^.\ ^y{-vi+.ii.+vi}-^^'(^«^^)- («^-i)(--^ 
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Putting these values in (9.45), multiplying through by (|^) and using 

dt 



(c.f. (9.32)), gives 

fdrA' d dfg 
\ dX J df dX 



dr 



a^ — l)(re — f) 1 (a^ — l)(re — f) 



8F^t F^Vi SF^t 

which yields 



±dTa ^ '\__ ( _ {a -l)(re-f) \ fdrA . . 

dfdX v^F^r 4v^ ]\dx)' ^""-^^^ 

and this implies (9.39) by easy estimates. Now (9.46) is a first order ODE 
for ^ with the initial condition ^(0) = 1. The initial condition fa(0) = 
then determines f^ from -jf. 



In the case a = 1, (9.46) reduces to 

d dfi 1 f dfi 



df dX \/t\ dX 



(9.47) 



which integrates exactly as follows. First set w = ^ so (9.47) becomes 



du 1 



df y/l 

then use v = y/t = |f + ^/t'e^ dv = ^df to get 

rdu r 2dv 



u /r V 



which, using Ug = 1, integrates to 



S-"4- f^-^^) 



Solving (9.48) gives t = (^f + v^) , dt = \ltdf giving 

ie,dX = yidi. 



Using Ag = 0, this integrates to 



3/2 



A = -V^e 
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(9.49) 



which upon solving for t gives t = ^i(A) in agreement with (9.40). Using 
t = (^f + ^/t^) in (9.40) then gives (9.41), thereby also implying (9.42). 



To obtain (9.43), write (9.46) in the form 
d dfa 



^ {l-0{l)|a-l|CJ^*''° 



df dX ^Ji 

Following the argument above then leads to 

dfn t 



dX 



(9.50) 



u 



= ^ + 0(l)|a-l|e 



dX t 
which directly implies (9.43). To see this, write (9.50) as 



(9.5i: 



(9.52) 



df y/t 

where we have incorporated |a — 1| into the 0(1). Integrating (9.52) gives 

r du _ r 2dv 

which leads to 



;i + o(i)Ce), 



V 



u 



> r-v 0(l)dv 



But we can estimate 

r 0{l)dv 



'Vt~e 



V 



r dv 

<0(1) / — <0(l)ln 



V 



VTe' 



SO 



> rv 0{l)dv 



<e^(^)^^i"7fe = 



V 



0(l)Ce 



= 1 + 



Vi- Vi 



0(l)Ce 



and expanding leads to 



/./j.^ < 1 + ^^i^om + 0(1) f^li^om I c 



e? 



which gives (9.51) because ^/i — y^ = f/2 = 0{re)- Integrating (9.43) gives 
(9.44) by easy estimates. This completes the proof of Lemma 17. D 

By a small angle perturbation of the baseline null geodesic we mean a null 
geodesic solution {t,f,4>) of (9.27)-(9.29) defined for < f < rg, satisfying 
the initial conditions at A = 0, 

(t,f,0)L^„ = (te,O,e) (9.53) 



^(*'''^)|^-- s 



+ 0{^) * 



A=0 



d\ 



+ 0(e2), , (9.54) 



A=0 



where we set the perturbation parameter e equal to 0e- To estimate such 
geodesices at fixed a, we now define a consistent O(e^) asymptotic approxi- 
mation for the small angle perturbation of the baseline null geodesic. That 
is, we find an approximate solution involving a function 0(A) which we show 
satisfies the null geodesic equations to order O(e^) when e = 0e, such that 
0(A) = e0(A). We accomphsh this by defining an ansatz which satisfies (9.27) 
and (9.28) to O(e^), and such that setting the leading order 0(e) term in 
(9.29) equal to zero, gives an equation for 0(A) which when satisfied, reduces 
(9.29) to 0(e2) as well. 

Our goal is now to prove the following theorem, which is a refinement of 
Theorem 14. 



Theorem 15 Fix a and set e = 0e- Then the following asymptotic approxi- 
mation is valid for the geodesic equations (9.27)-(9.29) as e ^ 0: 

f(A) = f„(A) + eV(A) (9.55) 

t{f) = ta{X) + eH{f), (9.56) 

0(f) = e0(f), (9.57) 
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where f{X), i{f) and (f){f) are 0(1) as e — > 0. More precisely, the approximate 
solution f = fa{X), t = ta{X), (defined in (9.34), (9-35)), solve the t and f 
geodesic equations (9.27) and (9.28) to order 0{e^), and the (j) geodesic equa- 
tion (9.29) to order e. Moreover, the e-order equation determined by (9.29) 
then determines (j) to third order in ( by 

Kre) = l-^3^Ce+0(l)|a-l|Cf as r, ^ 0, (9.58) 



Ce = 



re 



Said differently, if (j) solves the (jy-geodesic equation (9.29) to order e^^ then it 
must satisfy (9.58) in the limit e — > 0, thereby verifying that (9.58) is correct 
asymptotically in the limit e — > 0. 

Since e = 0e and 0(re) = 0o/0e, it follows from (9.4) that 

Ca = lim— = 0(re), 

and so (9.58) verifies the asymptotic relation (9.5). Thus Theorem 15 implies 
Theorem 14, and is more accurately a refined restatement of it. 

In the proof of Theorem 15 we will need the following lemma. 

Lemma 18 Assume < f < r^, te < t < to, (e ^ -^- Then the ansatz 
(9.55)-(9.57) implies the following asymptotic estimates: 

^(r)=^e|l + ^| +0{e^)=t,{l + 0{Q] as re,e^O, (9.59) 



— - = v^+0(e^) as re,e^O, (9.60) 

\Aji 



Proof: From (9.56), 



t=Qr + v^j +0{e')=t,[l + ~^ +0(e2), (9.61) 
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which imphes (9.59) for te <t < to. Differentiating (9.59) gives, 

= Vt + 0{e^) as VeC^O 
as claimed in (9.60). D 



It remains only to give the 

Proof of Theorem 15: To analyze small angle perturbations of the baseline 
null geodesies we need expressions for the Christoffel symbols (9.30). Since 
we are only looking to solve the geodesic equations (9.27)-(9.29) to order 
e, we only need values for (9.30) up to order e, but since the F's involve 
derivatives with respect to 0, this requires that we use values for gij in (9.26) 
up to order e^. On the other hand, the inverse matrix g^^ = g^^ enters 
(9.30) undifferentiated, so we need only compute g^^ to order e . Using that 
= 0e = 0(l)e, we can write (9.26) in the matrix form. 



/ -Ff 



9^J 







A p2^ 0^-1^212 a^-l, 



r 

a 4 e 







r — rpjrr^ 



\ 



\ 







a^-l, 



r — Vpjrr. 



tf2 + 0(e2) j 



+ 0(6^), (9.62) 



which to order e reduces to 



^ a 



9ii 







V 











a^-1, 



Fit 



r — rp]rr. 



a^-l, 







r — TpWr^ 



\ 



tf' 



+ 0(6^ 



(9.63) 



/ 



To compute the inverse g'^^ to order e, it suffices to take the inverse of (9.63), 
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and a straightforward calculation gives 






9 



^3 — 



F2 











V 







1 


a^-1 {re-r)re 

4^2 ^2 f. 


At^F^ r ^ 


1 



\ 



^Oie' 



(9.64) 



/ 



The next lemma records the relevant values of F^- to within errors necessary 
for evaluation of the 0(e) part of (9.29) in the limit as rg, e — > 0. 



Lemma 19 As rg, e — > we have 



p2 _ 

J- 11 — 


a2 - 1 re ^ Id 


2 tf'*' 2f'>-d4, 


r2 _ 

^ 00 — 


2f2ta.^(^«) + -^»' 


p2 _ 

J-01 — 


-^01' 


p2 _ 

■■- 02 — 


2>«(^^)- 


p2 _ 

J- 12 — 


^ + 0{e^), 



K)^ + ^n, 



(9.65) 

(9.66) 

(9.67) 
(9.68) 

(9.69) 



where the errors Ef^ are given by 



E 



11 



^00 



Vr 
0(l)\-^-J^e^O(e'] 



K = Oil] 



\a — l\rl 



6 + 0(6^ 



Proof: Let the second term in (9.30) involving derivatives of the metric be 
denoted 



{ij, k} = {-gij,k + gktj + gjk,i} 
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(9.70) 



To estimate these brackets we use the following easily verified identities: 
9'' = ^^^K^^ + 0{e% 



1 



r = 4+o(.^), 






A straightforward calculation using these together with gij given up to order 
e^ in (9.62), gives the following values for (9.70): 

{11, 2} = {a' - l)fr,0 - f^F'j {11, 1} = §{F'Jt 

{00,2} = |(F,2) {00,1} = |(F„)2 

{01, 2} = {01,1} = 1 (9-71) 

{12,2} = 2tf {12, 1} = 0(^) 

{02,2} = f2 {02, 1} = 0. 

For example, the caclulation of {11,2} above entails an interesting cancella- 
tion at the leading order as follows. By (9.70), 

{11, 2} = {-^11,2 + 2^21,1}, (9.72) 

and by (9.62), 

-gn,2 = 2"^r'J-^iF!)t, (9.73) 

4 d(p 

and 

^2 _ 1 

= -2^-—rl(^+{a^-l)rer(^. (9.74) 
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Adding (9.73) and (9.74) verifies {11,2} in (9.71). (Note the cancellation at 
leading order.) 

Using the values in (9.71) we can compute: 



where 

^11 - -^^0 71 17I\^a)9 + <^(e j, 



8F2 tr df 



from which (9.65) and (9.70) follow. 



where 

a^ — 1 (re — r)re d 
8F^ t^r df 

from which (9.66) and (9.70) follow. 



r^o = ^^{00, 2} + ^9^1 {00,1} 

^ '^■F!) + E-i„ (9.76) 



Ei, = \J'\'" iziK}'P + 0(<i% 



where 



r^i = ^^''{01, 2} + ^^^HOl, 1} = ^o\ (9.77) 



2 _ a^-l(re-f)re - .x 

^01 - ^;^^^ :^ + ^(e ), 



from which (9.67) and (9.70) follow. 



Tl = Ifi2^{12,2} + ^5^1{12,1} 

= l + 0{e') (9.78) 
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giving (9.69). And finally, 

verifies (9.68) and the proof of the Lemma 19 is complete. D 

To prove Theorem 15 we must show that (9.55)- (9.58) solve the geodesic 
equations (9.28)-(9.29) to order 0{e^). To start we verify (9.29). Using that 
(/) = 0(e), we can write (9.29) as 

_(f(i>_-n2 { dt\'^ ,T^2 {df\^ , op2 f dt\ fdf\ 

dx-' ~ ^00 [dxJ ^^iiKdx) '^ ^^oi[dx) \dx) 

+2r§2(f)(i) + 2r?,(S)(i) + o(^-^). 

the only derivative not appearing on the right hand side being ( ^ ) , which 
is O(e^). Multiplying through by (^) gives 

-(f)"'lS = rgo(f)' + r?, + 2rg,(|) (9.80) 

+2rg2(f)(f)+2r?,(f)+o(e^). 

Now since by (9.60) we know § = v^ + 0{e^), putting (9.65)-(9.69) with 
(9.70)-(9.70) in (9.80), using (9.60) gives 

\dx) dfdX~2Pd4>^^'^^^^ ^tf ^ 2Pd^^^^ 

+El + 2ElVi + 2 (^ + i) I + 0{e'). 

(Note the cancellation between the Pf^ and Fqo terms ^-^{F^).) 

Neglecting O(e^) terms, then multiplying through by \ftft^^ and using (9.70)- 
(9.70), gives the leading order equation for 

-Vir (S)-' 11 = ^^ + Vi (^ + 2) I + 0{l)|a - 1| {.} . 
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where 



-/. 






{•}=|^ + ^ + 7|=0(l)Ce+0(^')' 



(9.81) 



and we have used Lemma 18. Thus the leading order in e equation for is 



-r 



^/^(S) 



d d(f) 
df d\ 



^0 + v^(^ + 2)f + 0(l)|a-l|C 



(9.82) 



Now 



df\ d 
dX J df dX 



df\ d \ df 
dX I df \ dX df 



' df^ 
dX 



d ( df 



+ 



d' 



df \ dX \ df df^ ' 



and by (9.55), (9.39), 

f df\ d ( df\ 1 



{l + 0(l)|a-l|Ce} 



\dXj df [dX) ^/i 

Putting this into (9.82), neglecting higher order terms, gives 

-d^ 



^ {1 +««!«- i|«i- '-^rf,- 



-^^^-^(-A).- 



(9.83) 
+ 0{l)\a-l\Cl 



Now using that t as a function of f agrees with the baseline geodesic relation 

2 



1 f 



2.n 



t = tAi + 



that f I \lie is 0(l)Ce, and incorporating terms second order in Ce into 0(l)Cg, 
(9.83) is equivalent to 



, ^l± _ (g^-l) 
^^t, df^ " 



2 Ce0 + 2v^(l + i;^Hf + O(l)|a-l|e (9.84) 



Finally, we determine the constants Ci and C2 such that 



T IT 

[r) = 1 + Ci— + C2 -= 



+ o(i: 



r 



(9.85) 
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(With this ansatz, as well as equation (9.84) both reduce to functions of Ce 
at f = Tg.) Substituting this together with 




/ r 

+ 0(1] : 



e. 



into (9.84), evaluating at f = fg, and then collecting powers of Q, yields 

{2Ci} + jecs + "^ + Cij Ce + 0(1)4^ = 0. 

It follows then that (9.85) meets the 0-geodesic equation (9.29) to order 
0(l)Cg if we choose 

Ci = 0, C2 = -"^. (9.86) 

We conclude that the function 0(f) = e0(f) that solves the 0-geodesic equa- 
tion (9.29) to order O(e^), satisfies the asymptotic relation in Ce given by 

= i-^Cl + o(Cl). 

Since 0=1 solves (9.29) identically in the FRW case a = 1, (in the Stan- 
dard Model, the angle is constant along radial geodesies emanating from any 
center), it follows by smoothness that when (9.86) is assumed, the 0(1) in 
(9.85) is really 0(l)|a — 1|, and we have proven what was claimed in (9.58). 

To complete the proof of Lemma 15, it remains only to show that the ansatz 
(9.55)- (9.57) meets the t- and f-geodesic equations (9.27) and (9.28) to order 
e^ under the condition that (f) solves the e-order equation from the 0-geodesic 
equation (9.29) that we just analyzed. That is, it suffices to prove the follow- 
ing lemma: 

Lemma 20 Let ta{X) and ra(A) be as defined in the baseline formula (c.f. 
(9.34), (9.35)), and assume that the approximate geodesic 

(ta(A),f,(A),e0„(A)), (9.87) 

solves the (jy-geodesic equation to order O(e^). Then (9.87) also solves the t- 
and f-geodesic equations to order 0{e^). 

97 



Proof: Define the metric g, 

ds^ = -Fldt^ + Fltdf'^ = g.jdx'dx^, (9.88) 

so tliat ta(A), To (A) exactly solve the geodesic equations 

T-^O ^'^ HiJb , r^r\\ 



where 



Fl = F„^(C) = 1 + ^C^ (9.91) 



x2 



? = V' (9-92) 



r = Ve — r. (9.93) 

The point is that ( = r/Vi is obtained from ( = r jyji by setting cos = in 

the formula r = y r1 — 2ref cos0 + f^, thereby incurring an error no greater 
than O(e^). A direct consequence of (9.89)-(9.93) is that, 

Pi 

(9.94) 
(9.95) 

(9.96) 

Now to prove that (9.87) satisfies the equations (9.27) and (9.28) to O(e^), it 
suffices only to show that 

rg = f g. + 0(6^), z 7^ 2, jV 2, « = 0, 1 (9.97) 

V% = f^, + 0(e), J = 0,1, 2, a = 0,1. (9.98) 

(We don't need a = 2 because that applies only to the 0-geodesic equation, 
and we only need an order 0(e) on any term that gets multiplied by ^ 
because and its derivatives are 0(e) asymptotically.) 
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^ 2 

ur 


= 2r~ + 0(e2), 


dr 


= 0(e), 



Using (9.94)- (9.96) it is straightforward to obtain the following two tables 
giving the values of the brackets (9.70) associated with the full metric gij 
taken to order 0{e^) in (9.62). ({02,2} 



{00,0} = -^ {00,1} = ^ 

{11, 0} = -I {F^t} + 0(e2) {11, 1} = |(F^) + 0(e2; 

{01,0} = -^ {OlA} = UF'at) + 0{e'] 

{02,0} = -^ = 0(e) {02,1} = 0(e) 

{12,0} = 0(e) {12,1} = 0(e); 

{22,0} = 0(1) {22,1} = 0(1) 



(9.99) 



and the nonzero brackets of form {ij, 2} that are lower order than 0(e) have 
the following orders as e — > : 



{00,2} = 0(e) 


(9.100) 


{11,2} = 0(e) 


(9.101) 


{12,2} = 0(1). 


(9.102) 


Now it is also straightforward from (9.94)-(9.96) that 




Fi = F! + 0{^), 




'S^"S*°"''- 




?^f-'.'.. 




5 ^ o„. 





as e — )► 0, so we can replace F^ by F^ in (9.99), incurring no change in the 
stated errors in e. 

To verify (9.97), (9.98), consider first r^^-. Then 

r?, = lr{^J,0} = -If!{ij,0} = -If!{ij,0} + 0(e2). 

Thus it suffices to verify (9.97), (9.98) for {ij, 0} in place of T^j. But replacing 
F^ by F^ on the left hand side of (9.99) incurs no change in the e order errors, 
and this verifies (9.97), (9.98) in this case. 
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Consider next V}-. In this case 
'J 



rl; = ^ff"{u.i} + Jff'My.2}, 



where by (9.64). 



9'' = ^ = 0(1)' (9-103) 

9'' = -^^^^^^^4 = O {e) . (9.104) 



f'F^ 



r 



Therefore the result (9.97), (9.98) apphes after replacing F^ by F^ on the 
right hand side of (9.99), incurring no change in the stated errors in e. 

In short, the argument demonstrating that the t- and f-geodesic equations 
(9.27), (9.28) agree with (9.89), (9.90) to order e^ can be summarized as 
follows. The main observation is that by (9.63), the metric g in (9.88) agrees 
with g and its inverse g^^ agrees with g^^ to order e^, except in the gi2 
and g^'^ components, respectively, which both contain an order e error. As a 
consequence we need only show that the derivatives gij^k that differ by order 
e or by 0(1) from the corresponding derivatives in gij^k end up multiplied 
by terms 0(e) and O(e^) in the t- and f-geodesic equations (9.27), (9.28), 
respectively. 

Now the derivatives gij^k that differ by order e from the corresponding deriva- 
tives in gij^k are the t and f derivatives of ^'12, together with the derivative 
of a metric entry other than gu; and the only derivative gij^k that differs by 
0(1) from gij^k is gi2,2- Now derivatives of metric entries with respect to (f) 
can only appear as terms in F^- or Tj- with z = 2 or j = 2, or else as terms 
of the form g^^{ij,2} in Tjj, with i,j=0 or 1. But rj^2 is always multiplied 

by ^ in the geodesic equations, so these epsilon order errors become O(e^) 
in the equation, and can be neglected. Thus it remains only to see that the 
e-order derivatives p'12,0, 912,1 end up multiplied by order e in the t- and f- 
geodesic equations (9.27), (9.28), and that the zero order term ^'12, 2 ends up 
multiplied by order e^. But the derivative ^'12,0 appears only in r52 and rQ2, 
and gi2^i appears only in rj2 and P}^, and in each instance these derivatives 
are multiplied by either the order epsilon term g^'^ or else by the order 0(e) 
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term ^ in the t- and f-geodesic equations (9.27), (9.28); and the order one 
term ^'12,2 only appears in r22, which ends up multiphed by the O(e^) term 
as claimed. 






This completes the proof of Lemma 20, and the proof of Lemma 15 is com- 
plete. 



10 Concluding Remarks 

We have constructed a one parameter family of general relativistic expansion 
waves which at a single parameter value, reduces to what in this paper we 
call the FRW spacetime, the Standard Model of Cosmology during the radi- 
ation epoch. The discovery of this family is made possible by a remarkable 
coordinate transformation that maps the FRW metric in standard co-moving 
coordinates, over to Standard Schwarzschild coordinates (SSC) in such a way 
that all quantities depend only on the single self-similar variable ^ = f/i. 
Note that it is not evident from the FRW metric in standard co-moving co- 
ordinates that self-similar variables even exist, and if they do exist, by what 
ansatz one should extend the metric in those variables to obtain nearby self- 
similar solutions that solve the Einstein equations exactly. The main point 
is that our coordinate mapping to SSC form, explicitly identifies the self- 
similar variables as well as the metric ansatz that together accomplish such 
an extension of the metric. 

The self-similarity of the FRW metric in SSC suggested the existence of a 
reduction of the SSC Einstein equations to a new set of ODE's in ^. Deriving 
this system from first principles then establishes that the FRW spacetime 
does indeed extend to a three parameter family of expanding wave solutions of 
the Einstein equations. This three parameter family reduces to an (implicitly 
defined) one parameter family by removing a scaling invariance and imposing 
regularity at the center. The remaining parameter a changes the expansion 
rate of the spacetimes in the family, and thus we call it the acceleration 
parameter. Transforming back to (approximate) co-moving coordinates, the 
resulting one parameter family of metrics is amenable to the calculation of a 
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redshift vs luminosity relation, to third order in the redshift factor z, leading 
to the relation (1.1). It follows by continuity that the leading order part of an 
anomalous correction to the redshift vs luminosity relation of the Standard 
Model observed after the radiation phase, can be accounted for by suitable 
adjustment of parameter a. 

These results suggest an interpretation that we might call a Conservation Law 
Scenario of the Big Bang. That is, it is well known that highly interactive os- 
cillatory solutions of conservation laws decay in time to non-interacting waves, 
(shock waves and expansion waves) , by the mechanisms of wave interactions 
and shock wave dissipation. The subtle point is that even though dissipation 
terms are neglected in the formulation of the equations, there is a canonical 
dissipation and consequent loss of information due to the nonlinearities, and 
this can be modeled by shock wave interactions that drive solutions to non- 
interacting wave patterns. (This viewpoint is well expressed in the celebrated 
works [11, 6, 7]). Since the one fact most certain about the Standard Model 
is that our universe arose from an earlier hot dense epoch in which all sources 
of energy were in the form of radiation, and since it is approximately uniform 
on the largest scale but highly oscillatory on smaller scales^^, one might rea- 
sonably conjecture that decay to a non-interacting expanding wave occurred 
during the radiation phase of the Standard Model, via the highly nonlin- 
ear evolution driven by the large sound speed, and correspondingly large 
modulus of Genuine Nonlinearity^'^^ present when p = pc^ /?>^ c.f. [12]. Our 
analysis has shown that FRW is just one point in a family of non-interacting, 
self-similar expanding waves, and as a result we conclude that some further 
explanation is required as to why, on some length scale, decay during the 
radiation phase of the Standard Model would not proceed to a member of 
the family satisfying a 7^ 1. If decay to a 7^ 1 did occur, then the galaxies that 
formed from matter at the end of the radiation phase, (some 379, 000 years 
after the Big Bang), would be displaced from their anticipated positions in 
the Standard Model at present time, and this displacement would lead to a 
modification of the observed redshift vs luminosity relation. In short, the 



^^In the Standard Model, the universe is approximated by uniform density on a scale of a billion light years 
or so, about a tenth of the radius of the visible universe, [18]. The stars, galaxies and clusters of galaxies 
are then evidence of large oscillations on smaller scales. 

^■^ Again, Genuine Nonlinearity is in the sense of Lax, a measure of the magnitude of nonlinear compression 
that drives decay, c.f., [11]. 



102 



displacement of the fluid particles, (i.e., the displacement of the co-moving 
frames in the radiation field), by the wave during the radiaton epoch leads to 
a displacement of the galaxies at a later time. In principle such a mechanism 
could account for the anomalous acceleration of the galaxies as observed in 
the supernova data. Of course, if a 7^ 1, then the spacetime within the ex- 
panding wave has a center, and this would violate the so-called Copernican 
Principle, a simplifying assumption generally accepted in cosmology, at least 
on the scale of the wave (c.f. the discussions in [17] and [1]). Moreover, if our 
Milky Way galaxy did not lie within some threshold of the center of expan- 
sion, the expanding wave theory would imply unobserved angular variations 
in the expansion rate. In fact, all of these observational issues have already 
been discussed recently in [2, 1, 3], (and references therein), which explore 
the possibility that the anomalous acceleration of the galaxies might be due 
to a local void or under-density of galaxies in the vicinity of the Milky Way.^^ 
Our proposal then, is that the one parameter family of general relativistic 
self-similar expansion waves derived here are possible end-states that could 
result after dissipation by wave interactions during the radiation phase of the 
Standard Model is completed, and such waves could thereby account for the 



^'*The size of the center, consistent with the angular dependence that has been observed in the actual 
supernova and microwave data, has been estimated to be about 15 megaparsecs, approximately the distance 
between clusters of galaxies, roughly 1/200 the distance across the visible universe, c.f. [1, 2, 3]. 
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appearance of a local under-density of galaxies at a later time. 
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In any case, the expanding wave theory is testable. For a first test, we propose 
next to evolve the quadratic and cubic corrections to redshift vs luminosity 
recorded here in relation (1.1), valid at the end of the radiation phase, up 
through the j9 « stage to present time in the Standard Model, to obtain 
the present time values of the quadratic and cubic corrections to redshift vs 
luminisity implied by the expanding waves, as a function of the acceleration 
paramter a. Once accomplished, we can look for a best fit value of a via com- 
parison of the quadratic correction at present time to the quadratic correction 
observed in the supernova data, leaving the third order correction at present 
time as a prediction of the theory. That is, in principle, the predicted third 
order correction term could be used to distinguish the expanding wave theory 
from other theories (such as dark energy) by the degree to which they match 
an accurate plot of redshift vs luminosity from the supernove data, (a topic 
of the authors' current research). The idea that the anomalous acceleration 
might be accounted for by a local under-density in a neighborhood of our 
galaxy was expounded in the recent papers [2, 3]. Our results here might 
then give an accounting for the source of such an under-density. 



^^Thc following back of the envelope calculation from [16] provides a ballpark estimate for what we might 
expect the extent of the remnants of one of these expanding waves might be today if our thesis is correct 
that wave interactions and dissipation by strong nonlinearities during the radiation phase were the primary 
mechanisms involved in formation of the wave. For this, note that matter becomes transparent with radiation 
at about 300,000 years after the Big Bang, so we might estimate that the wave should have emerged by about 
tendrad ~ 10^ ycars after the Big Bang. At this time, the distance of light-travel since the Big Bang is about 
10"'' lightyears. Since the sound speed c/v3 « .58c during the radiation phase is comparable to the speed 
of light, we could estimate that dissipation that drives decay to the expanding wave might reasonably be 
operating over a scale of 10^ lightyears by the end of the radiation phase. Now in the p — expansion that 
follows the radiation phase, the scale factor (that gives the expansion rate) evolves like 

[18], so a distance of 10'"' lightyears at t = tendrad years will expand to a length L at present time tpresent ~ 
lO^'^ years by a factor of 

^KJ-present) ^^ \ ) 

R{tendrad) (10^)^^^ 

It follows then that we might expect the scale of the wave at present time to extend over a distance of about 

i = 5 X 10^ X lO'^ = 5 X 10^ lightyears. 

This is a third to a fifth of the distance across the visible universe, and agrees well with the extent of the 
under-density void region quoted in the Clifton-Ferriera paper, with room to spare. 
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= lO''-^ > 5 X W. 



In summary, the expanding wave theory could in principle give an explana- 
tion for the observed anomalous acceleration of the galaxies within classical 
general relativity, with classical sources. In the expanding wave theory, the 
so-called anomalous acceleration is not an acceleration at all, but is a correc- 
tion to the Standard Model due to the fact that we are looking outward into 
an expansion wave. The one parameter family of non-interacting, self-similar, 
general relativistic expansion waves derived here, are all possible end-states 
that could result by wave interaction and dissipation due to nonlinearities 
back when the universe was filled with pure radiation sources. And when 
a 7^ 1 they introduce an anomalous acceleration into the Standard Model of 
cosmology. Unlike the theory of Dark Energy, this provides a possible expla- 
nation for the anomalous acceleration of the galaxies that is not ad hoc in the 
sense that it is derivable exactly from physical principles and a mathemat- 
ically rigorous theory of general relativistic expansion waves. In particular, 
this explanation does not require the ad hoc assumption of a universe filled 
with an as yet unobserved form of energy with anti-gravitational properties, 
(the standard physical interpretation of the cosmological constant), in order 
to fit the data. The idea that the anomalous acceleration might be accounted 
for by a local under-density in a neighborhood of our galaxy was expounded 
in the recent papers [2, 3]. Our results here might then give an accounting 
for the source of such an under-density. 

In conclusion, these new general relativistic expanding waves provide a new 
paradigm to test against the Standard Model. Even if they do not in the end 
explain the anomalous acceleration of the galaxies, one has to believe they 
are present and propagating on some scale, and their presence represents an 
instability in the Standard Model in the sense that an explanation is required 
as to why small scale oscillations have to settle down to large scale a = 1 
expansions instead of a 7^ 1 expansions, (either locally or globally), during 
the radiation phase of the Big Bang. 
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